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1. INTRODUCTION 

In attempting to formulate a given "equilibrium-type" of continuous 
system in mechanics one may either use some kind of a variational princi- 
ple and reduce the problem to a minimization problem subject to certain 
constraints, or, as is more commonly the case, one may directly apply the 
equilibrium principles and reduce the problem to a boundary value problem 
which consists of a (system of) differential equation(s) subject to cer- 
tain boundary conditions. Even though in most practical applications the 
minimization problem is further reduced to a boundary value problem, it 
may also be solved directly by using an approximate technique such as the 
Ritz's method. To facilitate the definition of certain concepts consider 
the following boundary value problem in two dimensions: 

^ 2 ^( 0 ) = f(xj,X2> , (xj, X2 )eD , (1.1) 

B^(u) = §i(s) , (i = l,..,m) , seS (1.2) 


where is a differential operator of order 2m, Xj^,X 2 are the spatial 
coordinates, D is the domain of definition for the unknown function u, S 
is the boundary of D, B^ (i = 1, ..,m) is a differential operator (contain- 
ing u and its normal derivatives) of order (at most) 2m- 1, f and g are 
known functions, and s is any convenient coordinate defining the point on 
the boundary (say, the arc length). The domain D may contain the point at 
infinity and may be multiplyconnected. Contours forming the boundary are 
assumed to consist of piecewise smooth arcs. The points on S at which the 
tangent has a discontinuous slope will be called the points of geometric 
singularity . 
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There is another type of singular point on the boundary which results 
from the change in the nature of the homogeneous operators B^(u) specify- 
ing the boundary conditions. Such a point on a smooth boundary either 
side of which at least one of the operators (i = l,,.,m) has different 
behavior is called a point of flux singularity . Note that if the behavior 
of an operator B^ changes at a "corner point" of the boundary, this point 
is then a point of both a geometric as well as a flux singularity. If 
the homogeneous operators B^ remain unchanged on each closed contour (but 
not necessarily the same on all contours) , then the corresponding problem 
is called an ordinary boundary value problem . On the other hand if there 
are points of flux singularity on the boundary, the problem is called a 
mixed boundary value problem . 

In working with mixed boundary value problems it is often advantag- 
eous to keep in mind that the physical system has generally two types of 
quantities, namely, the potential and the flux type quantities. In poten- 
tial theory the meaning of these concepts is unambigious and clear. They 
are, for example, identified by temperature, velocity potential, electro- 
static potential, mass concentration, or displacement (in anti-plane shear 
problems) as the potential type quantities, and heat flux, velocity, elec- 
trostatic charge, mass rate of diffusion, or stress, as the corresponding 
flux type quantities. Similarly, in solid mechanics one may classify the 
displacements and the stresses (or the strains) as respectively the poten- 
tial and flux type quantities. The physics of the problem requires that 
the "potential" be bounded and continuous everywhere in IH-S> including 
the points of both geometric and flux singularity. 

To fix the ideas, consider the following simple problem in a wedge- 
shaped domain: 

1 


V^u(r,9) = 0 , 


0 <r<«> , 


o<e<0 


(1.3) 


where the known functions f and £2 are such that the global equilibrium 
is satisfied and 0<6^<2 tt. The problem is an ordinary boundary value 
problem and r = 0 is a point of geometric singularity. If in the neigh- 
borhood of r = 0 fj^ and £2 are zero, simple application of Mellin trans- 
form would indicate that for small values of r the components of the flux 
vector are of the following form: 


f, 

■ "I 1 

U I 


- = Kr"” “o^' “o sin 7 ^^ + 0 (r^^'^"®o^^®o) , 


J, 3u _ ^ (iT-e ^)/0 
r 30 


Jl9 

9 


o 

3u ^ (7T-0 )/0 7r0 , (27T-0 o)/0on 

TT- = -Kr o' ' o cos - 5 —+ 0(r^ °) , 

dr Oq 


(1.5) 


where K is a constant. Note that at r = 0 the flxix becomes unbounded for 

7r-0 

7r<0^<27T, the corresponding power of singularity being 0 > ’ a ~~ ° >- 0 . 5 . 

°o ~ 

If the boundary conditions (1.4) are replaced by 


u(r, 0 ) = 0 , 




0<r<jo 


( 1 . 6 ) 


it is seen that the problem is a mixed boundary value problem in which 
r = 0 is a point of both geometric and flux singularity. Hence, in this 
problem one would expect a stronger flux singularity than in the previous 
problem. The asymptotic solution of the problem for small r may again 
be expressed as 


u(r, 0 ) - Kr^^^^o sin 




1 du _ „ 7T ^.(Tr-29^)/29 
r 30 ^ 20„ ^ 


o cos 


7T9 

26„ 


+ O(r«’'-2eo)/20o) , 
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(1.7) 


= If _!!_ ^('n‘-20o)/20o 

9r ■ 20 ^ 

o 


sin 


TT0 

20 

o 


+ 20q)/20q^ 



From (1.7) it is seen that for 'n'/2<0 <27T the corresponding power of the 

7T-20O 

flux singularity is 0> — — ^-3/4, which is stronger than that found in 

o 

the previous problem for the same angle 0^. 

Let us now consider a special case of this problem in which 0 q=7T. 

Here the boundary is the infinite line; consequently the geometric singu- 
larity is removed but r=0 remains to be a point of flux singularity 
having a power —1/2 . As will be seen later in this article, - 1/2 power • 
singularity is quite typical for the points of flux singularity on a 
smooth boundary. It will also be seen that however, somewhat contrary 
to the general expectation, this is not always the case in mixed boundary 
value problems With boundary conditions containing the potential and the 
flux. 

Finally, consider the following (mixed) boundary conditions 

u(r,0) = 0 , a<r<b , 

u(r,0) = g(r) , 0<r<a , b<r<“ , 

u(r,0Q) = f (r) , 0<r«» . (1.8) 

where f and g are knoxra functions. It is seen that r = 0 is a point of 
geometric singularity and the points r -a and r = b on the smooth boundary 
(0 = 0, G<r<“) are points of flux singularity. The problem is a mixed 
boundary value problem. Around the geometric singularity r =0 the solu- 
tion is expected to behave as in (1.5) and around the points of flux singu- 
larity it will have a behavior as in (1.7) X'J'ith = F ^’^'1 necessary coordi- 
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nate transformations. It should again be emphasized that the standard 
-1/2 power of the singularity at r = a>0 will be a function of 6^ for 
a = 0 which, depending on the value of 9 q, may be stronger or weaker 
than - 1/2 . A similar phenomenon will be discussed in connection with 
a contact problem in elasticity later in this paper. 

2, DEFINITIONS: MULTIPLE SERIES EQUATIONS , MULTIPLE INTEGRAL EQUATIONS 

In considering the solution of a given mixed boundary value problem 
perhaps the simplest technique is the direct application of the method 
of complex potentials provided the problem admits such potentials and 
the domain and the boundary conditions are suitable for such an applica- 
tion. In this case the problem is reduced to a Riemann-Hilbert problem 
for a (system of) sectionally holomorphic function(s) which may be 
solved in a straightforward manner. On the other hand if one applied a 
more standard technique such as, the separation of variables, integral 
transforms, or the method of Green's function, the mixed boundary condi- 
tions invariably lead to a formulation involving "dual series equations", 
"dual integral equations" or "singular Integral equations". Again, to 
facilitate the basic understanding of these notions, the definitions will 
be preceded by the formulation of some simple examples. 

2.1 Multiple Series Equations. 

As a first example consider the following mixed boundary value ptob- 
lem in potential theory for the unit circle: 

V^u(r,6) = 0 , 0<r<l , 0<9<2 t , (2.1) 

u(l,9) = g(9) , 9eL^ , (2.2) 

u(l,9) = f(9) 5^^ : (2.3) 


where 


N 

= 2 = (r=l; a^<e<b^) , (2.4) 

and is the complement of on the unit circle. Using the technique 
of the separation of variables the solution of (2.1) may be expressed as 

u(r, ) = + E r (A^ cosn6 + sin n 0) . (2.5) 

Formally, substituting from (2.5) into the boundary conditions one obtains 
the following system of equations to determine the coefficients A^ and B^: 

CO 

E n(A^ cos n0 + B^ sinn0) = g(0) , 0eLj^ , 

CO 

A + E (A cos n0 + B sinn0) = f(0) , 9 eL„ . (2.6a,b) 

o^^n n / 2 

In the case of Dirichlet (L^ = 0) or Neumann (L 2 = 0) problems (2.6b) or 
(2.6a) give the unknown coefficients directly by expanding f(0) or g(6) 
into Fourier series in (0,2^). However, in the present problem the func- 
tions sinn0 and cos n0 (n = 0, 1,2, . . . ) are not orthogonal on and L 2 
and hence, (2.6) is at best equivalent to (or can be transformed into) 
an infinite system of algebraic equations. 

The structure of (2.6) is typical of the mixed boundary value prob- 
lems defined in a bounded domain a<x<b which may be expressed as 


00 


E A k, (n,x) 

1 T1 i 

= fj^(x) , 

xeL^ , 


00 

E A k„(n,x) 
, n 

= f2(x) , 

x£L2 , 

(2.7) 


or, more generally, for J sets of coefficients A^^, (j - 1, . . ,J; n = 1,2, . . . 
one obtains 


6 


S ^ S k (n.x)A^ - (x) , 

j = l n=l J 


xeL, 


(i=l,..,J) , 


2 E (n,x)A = (x) , 

j=l n=l 


xeL 


2 ’ 


(i=l,..,J) , (2.8) 


where again 


N 

^ ^Ik ’ 


^Ik " 


(s-jh) , 


W“k+i • 


a<aj^ 


b«5b , 


(2.9) 


-1 i 

and the input functions f or f and the kernel functions k or k . are 

r r r rj 

known. The system of equations such as (2.6), (2.7) or (2.8) are defined 
as dual series equations . 

Going back to the problem for the unit circle (2.1), if one assumes 
that the boundary L = (0,27 t) is divided into three parts with the following 
boundary conditions'^ ^ 


-^u(l,9) = f^(0) , 


0£L, 


u(l, 0 ) = f 2 ( 0 ) 


0 eL, 


h^u(l, 0 ) + h 2 ^^(1.9) = f 3 ( 9 ) 


0eL. 


( 2 . 10 a-c) 


where L , (r= 1,2,3) is the union of nonintersecting arcs L (r = 1,2,3; 

X Xj 

j = 1,2,..,J^) on the unit circle with + L 2 + = L = (0,2tt). h^^ and 

h 2 may be functions of 9. Again, formally from (2.1) and (2.10) it follows 
that' 


E n(A cos n0 + B sin n0) = f , (9) , 

j n n 1 


9eL. , 


( 1 ) 


This; is a problem, for example, in heat conduction in which, in addition 


to specifying the heat flux and temperature on parts of the boundary and 
L 2 J there is free convection taking place along L 3 where h. is the coeffic- 
ient of heat convection, hg is the coefficient of heat conduction, and f^ 
is related to the environmental temperature u^ through hjiUco = f^- 
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A + 2 (A cos n0 + B sinn0) = f^C©) > 

o , n n 2' ' ’ 


0£L, 


h,A + h, E (A cos n0 + B sinnO) 
lo l^ n n 

00 

+ h 2 E n(A^ cos n 0 + B^ sinnO) = £ 3 ( 0 ) » * ( 2 . 11 a-c) 


In this problem as seen from (2.10) the boundary conditions are defined in 
terms of three distinct operators on three separate parts of the boundary 
giving rise to a set of three series equations described by (2.11). Thus, 
these equations will be called triple series equations . In general, then 
a set of series equations for a system of unknown coefficients Aj^, A 2 >.-. 
of the form 


E A k. (n,x) = f . (x) , 


M 

xeL . , E L . = L , 

X ’ 1 X 


(i - 1,..,M) , 


will be called multiple series equations (of multiplicity 


( 2 . 12 ) 


In the mixed boundary value problems described by (2.6) and (2.11) the 
intersections of the boundary segments are points of "flux singularity". 
The discussion given in the previous section for the wedge would indicate 
that at least at some of these points the flux will be unbounded, at others 
the behavior of the solution is not known beforehand. Hence, the technique 
developed to solve the multiple series equations will not only have to be 
sufficiently general to apply to a great diversity of problems but will 
also have to lend itself to the correct treatment of the singular nature 

^^^Note that this definition is different than that found in current liter- 
ature (e.g., [ 1 ]) on "dual series" and "dual integral" equations where the 
multiplicity of the equations is taken to be the number of independent seg- 
ments on the boundary rather than the number of independent operators de- 
fining the boundary conditions. 
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of the solution. Aside from the method of complex potentials whenever 
applicable, it appears that the method of singular integral equations 
is the only approach which fulfills these requirements. Because of the 
singularities, since the infinite series giving the components of the 
flux vector will be divergent at certain points on the boundary, it is 
clear that any direct method reducing the multiple series equ 'itions to 
an infinite system of algebraic equations in the unknown coefficients 
(which has no closed form solution) will not be acceptable. 

2.2 Multiple Integral Equations 

Consider now the equivalent problem in potential theory for the half 
plane y>0 . The problem is stated as follows : 


V^u(x,y) = 0 , 


u(x,0) = fj^(x) 


u(x,0) = ^2^^^ ’ 


(_oo<x<oo ^ 0<y<«>) , 


xeL, 


xeL. 


N , 

L, = 2 L-, . , L, . = (a.,b.) , a.<h.<a.,, , 

1 ^ li ’ li x’ X X X x+1 


-oo<a 


1 ’ 




(2.13) 

(2.14) 

(2.15) 


(2.16) 


where L 2 is the complement of on (-°°<x<‘») , fj^ and f 2 are known functions 

2 2 

and are such that u-»-0 as x + y (i.e., any homogeneous solution behaving 
differently at infinity has been separated) . Using Fourier transforms, 
the solution of (2.13) may be expressed as 


u(x,y) = / 


-y a -iax 
i ' da , 


(2.17) 


where A (a) is an unknown function. Substituting from (2.17) into (2.14) 
and (2.13) formally we find 


|a|A(a)e ^ ^ da = f ^^(x) , xeLj^ , 


A(a)e~'^“'' da = f2(x) , 


xcL, 


(2.18a,b) 


If or 1,2 is zero, (2.18) may be solved in closed form in terms of inver- 
sion integrals. Integral equations of the form (2.18) or more generally 
the pair of integral equations 


/ kj^(x,a)A(a)da = fj|^(x) , xeLj^ , 

ii 

/ k2(x,a)A(a)da = > xeL2 , 


L^+L2=L , 


(2.19) 


in which the kernels k^ and IC 2 are different, are called a set of dual 
integral equations for the unknown function A(a) . If the problem Involves 
more than one unknown function, the boundary conditions would lead to a 
system of dual integral equations for the unknown functions A^ (a) , (j=l,.,,J) 
of the following form; 

J 


/ E k^. (x,a)A.(a)da = f^(x) 
L 1=1 ^ ^ 


xcIj ^ If 


/ Z 14. (x,a)A (a) da = f^(x) 
L 1 ^ 3 


xeL^ 


i — 1, . . , J) 


+ L 2 = L . 


( 2 . 20 ) 


In this problem too one may consider the following more general boundary 
conditions: 


ay 


u(x,0) = fj^(x) , xeLj , 


u(x,0) = f2(x) , 


xeL 


2 ’ 


h^u(x, 0) ^2 ~ ^3^^^ ’ ^^^3 * 


(2.2la-c) 
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where Lj^ + L2 + = (-“,«>) , and L^, 

ing line segments L„ (k = 1,.,,K,) 

XIC ' X 

(2.21) it follows that 


(i = 1,2,3) consists of nonintersect- 
on the real line. From (2.17) and 


- |alA(a)e da = > xeLj^ , 

A(a)e da = f2(x) , XEL2 j 

CO . 

/ (h, -h„ |aj )A(a)e da = f«(x) , xeL^ r (2.22a-c) 

—00 1 /. O jj 


Equations (2.22a-c) form a set of triple Integral equations for the unknown 
function A (a) . More generally 


/ S k’^. (x,a)A. (a) da 
LI J 


= f^(x) 
m 


M 

Z L = L , n = i, 

-1 ™ 
m=l 


j 


,J 


xEL , 
m ’ 


y 


m 



(2.23) 


is called a system of multiple integral equations (of multiplicity M) for 
the unknown functions Aj^,..,Aj. 


As in multiple series equations, in problems formulated in terms of 
multiple integral equations the singular nature of the solution is gener- 
ally not known beforehand. Therefore, in these problems too it is import- 
ant that the method of solution developed to solve the integral equations 
be not only sufficiently general and effective but also give the correct 
behavior of existing singularities. In this respect, particularly in dealing 
with somewhat unusual niixed boundary value problems, the methods of complex 
potentials and singular integral equations appear to be far superior to 
the standard operational techniques. An extensive treatment of the opera" 
tional techniques for the solution of dual series and dual integral equa- 
tions may be found in a recent book by Sneddon [1]. [2-13] are some of 

the outstanding references on the theory and applications of the complex 
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potentials and the singular integral equations. In this article the 
primary emphasis will he on the recent developments concerning the 
methods of solution of the singular integral equations and particularly 
the applications to some mixed boundary value problems with uncommon 
singularities . 

3. APPLICATION OF COMPLEX POTENTIALS 

In this section the direct application of complex potentials will 
be described by considering some relatively simple examples. 


3.1 A Problem in Potential Theory 

Consider the mixed boundary value problem in potential theory for 
the half plane (-“<x <“5 y>0) which is formulated by (2. 13)-(2 . 16) (Figure 1) 
Let the harmonic function u(x,y) be the real part of a complex potential 
F(z), z = X + iy. F(z) is holomorphic in the upper half plane S where 
the derivatives of u may be expressed as 


2 1^ = F'(z) + F'(z) , - 2i 1^ = F'(z) - F’(z) . 


(3.1) 


Noting that if zeS and z->-t + iO then zeS and z-^-t-iO, taking the boundary 
values of (3.1) * and using the conditions (2.14) and (2.15) we obtain 


.+ 


F’ (t) - F' (t) ^ - 2if^(t) , 


t£L, 


F*+(t) + F’ (t) = 2f’(t) , t£L2 


(3.2a,b) 


where (2-. 15) is used in differentiated form which means that for single- 
valuedness the solution must satisfy the following conditions: 
b, 




(3.3) 
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where k = if a,>-°o, h„<°°, and k = 2,..,N-1 if a, = and b., = 

I N 1 N 

- 1 - — 

Since F'(z) is holomorphic in S , F'(z) will be holomorphic in S . 


If is finite we define a new sectionally holomorphic function by 
G(2> = i 


^F'(z) , zeS"^ 


V 


F’(z) , zeS 


(3.4) 


From (3.2) and (3.4) it then follows that 


G (t) + G (t) = - 2ifj^(t) , » 


G^(t) - G“(t) = 2f^(t) , t£L2 


(3.5a, b) 


From (3.4) and (3.5) it is seen that G(z) is holomorphic everjTwhere in 
the complex plane except on and on that part of L2 on which f^^t) is 
not zero. Depending on the behavior of u, G may also have a pole of 
finite degree at infinity. For example, if there is a uniform "flux" at 
infinity given by 


9u , . 9u , . 

+ ^ * Pi + 


3x 


9 2 

(X +y ->«>) , 


(3.6) 


we have 


lim F' (z) = p, - ip. 


\z\-^ 


(3.7) 


Consider now the related homogeneous Riemann-Hilbert problem given by 


X'^(t) + X"(t) =0 , tsL. 


X'‘'(t) - X (t) = 0 , teL, 


N 

^1 " j ^Ik 


(3.8) 


where X(z) is the fundamental solution of the original Riemann-Hilbert 
problem (3.5) which is clearly determinate within an arbitrary multiplica- 
tive analytic function. Referring to, for example, [3], the general solu- 
tion of (3.8) may be expressed as 
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N 1 ^ 1^0 

+ Civ ^ , vTT + 3r 


X(z) = P(z) fTI(z-aj^) 2 “k(z-b^)2 *^k 


(3.9) 


where and 3j^(k == 1,..,N) are arbitrary integers (positive, negative, 
or zero) and P(z) is an arbitrary polynomial. At this point it should 
be strongly emphasized that in (3.9) as well as in the applications of 
the function-theoretic methods to the singular integral equations else- 
where in this article the arbitrary integers oi^, 3^ cannot be determined 
from purely mathematical considerations. To do this the physics of the 
problem has to be properly taken into account. In the present problem 
the points a^ and b^ are the type of singular points at which flux vector 
has an integrable singularity. Therefore, otj^ = 0, 3^ = -1, (k = 1,..,N), 
and ignoring the arbitrary polynomial, the fundamental solution of the 
problem becomes 


N 


X(z) = T71 (z-a^) 2 (z-\) 2 


(3.10) 


where for the particular branch considered lim z^X(z) =1. 

j z|-^ 

+ 

Dividing both sides of (3.5) by X (t) and using (3.8) we find 


. G(t). + _ .G(t) .~ ^ 

^X(t)^ ^X(t)'^ 


r 




- 2if^(t)/x^(t) 

. 2f’(t)/x'^(t) , 

\ 2 


teL 


1 


teL, 


(3.11) 


Equation (3.11) is now a simple boundary value problem the general solu- 
tion of which having a finite degree at infinity may be written as [3], 


G(z) 

X(z) 


- 


2if^(t) 


^1 (t-z)x'‘'(t) 


dt + 




2 f^(t) 


27ti L,) / 

2 (t-z)X (t) 


dt + \(z) (3.12) 


where P^(z) is an arbitrary polynomial of degree It. Erom (3.4) and (3,7) 
it is seen that G(z) has a pole of order zero at Infinity. Thus, from 
(3,10) and (3,12) it follows that the degree of is N, and (3.12) becomes 
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S(,) = - ^ 


„+,.. TTi •'L 


N 


1 (t-z)X'(t) 


f’(t)dt I, , 

+ X(z)£c^z . (3.13) 




2 (t-z)X'(t) 


where are arbitrary constants. From (3.4), (3.7) and (3.13) it 


is easily seen that 


^1 ~ ^^2 


(3.14) 


Noting that is finite, the remaining N constants c^, . . are determined 

from the single-valuedness conditions (3.3). From 


2 ^ u(x,+0) = F''*'(x) + F'“(x) = g‘*‘(x) - G“(x) 


(3.15) 


these conditions may he expressed as 


j [G'^(t) - G~(t)]dt = f2(b,p - f2(V » (k = lv..,N) (3.16) 

k ■ 


giving, with (3.13), a system of N linear algebraic equations in c^,.. ,c^_j^. 
This completes the solution for finite L^. 


If l,^ instead of is finite (i.e., = “) the procedure to 


solve the problem is quite similar to that given above with the following 
main differences: the branch cut should be introduced along 1-2 defining 

G(z) as 


F'(z) , zef3 


+ 


G(z) - 


F’(z) 


(3.17) 


zeS 


which would give the fundamental solution as follows: 


N-1 


X(z) = |~1 (z-b,)"^/2(z-r 


(3.18) 


Considering only the solution for which 


00 3 

u(x, 0) dx = Q = finite 


(3.19) 
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and noting that X(z)-^z^ for Jz [->«>, it may be seen that in the expression 
of G(z) similar to (3.13), the degree of the arbitrary polynomial will be 
N-2, i.e., the solution will contain N-1 arbitrary constants c ,..,c„ 

These constants are then determined from (3.19) and (3.3), completing 
the solution. 

3.2 The Case of Periodic Guts 

In the problem considered in the previous section let the cuts (a, ,b.) 

iC xC 

be equal in length and be equally spaced. Thus, one may define the end 
points of 2m+l cuts by 

a^^ = 2kb - a , bj^ = 2kb + a , (k = 0,+l,..,+m) 

m 

^Ik " ’ ^1 "-m ^Ik • (3.20) 

The fundamental solution (3.10) of the problem then becomes 

X(z) = (z-aj^)“^^^z-bj^)"^/^ 

= A[(z+a) r^(l-(|g)^(z-a)rj(l-(|^^ (3.21) 

where A is a constant given by 
A = (-l)“(2b)"^™(m:)"^ 

Because of the homogeneous nature of (3,8), since X(z) is determinate within 
a multiplicative analytic function, the constant A may be (and will be) 
ignored. It is seen that if we now let nrH», geometrically the problem 
becomes that of a plane with periodic cuts. In addition to this if the 
functions f^^ and f2 are assumed to be periodic then we have a problem for 
a half plane with periodic mixed boundary conditions. For simplicity, here 
it will be assumed that any homogeneous "loading" condition at infinity has 


been separated through a proper superposition and consequently in the 
problem under consideration F(z)->0 as First ignoring the constant 

A and then using the relation [14] 

oo 2 

sin0 =0 ] I (1 2~2 ^ ’ (3.22) 

in limit from (3.21) X(z) may be obtained as 

■y / ^ -1 ^ , 2 TTa V "l/2 

X(z) = (sin - sxn , (3.23) 

where the particular branch which is positive for a<Re(z) = t<2b-a will 
be considered . 

With X(z) as given by (3.23) the solution (3.12) is still valid. 

Since the geometry and the boundary conditions in the problem are periodic 
in X with a period of 2b, the potential G(z) must also be periodic in z 
with a real period 2b. On the other hand from (3.23) it is seen that 

X’^(t) = - X"(t) = - X'‘'(t+2b) = X"^(t+4b) , teL^ , 

X'^(t) = X"(t) = - X'^(t+2b) = X'*'(t+4b) , teL^ (3.24) 

namely, X(z) is periodic with a period of 4b. Therefore, in order to have 
a periodic potential with (real) period 2b, the arbitrary polynomial in 
(3.12) must be of the following form 

Pj^(z) = (Bj^ sin + B 2 cos '|^)Q(z) (3.25) 

where Q(z) is an arbitrary periodic analytic function with period 2b and 
Bj^ and B 2 are arbitrary constants. From 

lira G(z) = 0 = lim X(z)P^(z) = (Bj+iB 2 )Q(“) ^ ^ ^ (3. 26) 

|z|-x=° , [z|^ 

' 17 ' . ■ ■ ' ■ 


it follows that Q(“) = 0, or since Q is analytic in the entire plane 
Q(z) - 0. 


The solution of the problem may then be written as 


I. A 


+ X(z) “ j.\+l 


k (t-z)X (t) 
where X(z) is given by (3.23). Letting 




t == T+2kb , fj^Ct) = fj^(t+2kb) = , 




f;(t) = goCT) ,, X (t) = X(t) , -b<T<b , k= 0,4-1,. 


and using (3.24) the integrals in (3.27) may be modified as follows; 


1 “-“-'"k (t-.)xV) -»(T-.)«kb 

gl(T)dT , 

-•a X(T) T-z 1 (T-..z)2_(2kb)^ 

In (3.29) the series may be summed by using [14] 


1 . “ 20 ^ ” 20 _ _ _ Q 

0 J .2 2 ^cotirO , E 2 2n-1.2 .2 ® 

1 0 -n 1 IT ( — 2 — ) 


which gives 


cr (x) dT 

h = ^ ^a^x(ir i ^ 


7t: ra 
2b -a 


gl('r)dT 


sin -^(t-z) 


Similarly 


” A+IJ12 

-w b 


f I ( t) dt 




k (t-z>}T (t) 


JL 

® X(T)sin -^(T-a) 


(3.27) 


(3.28) 


(3.29) 


(3.30) 


(3.31) 


(3.32) 






Thus, the solution becomes 


- X(z) j-a X(z) j2b-a 

X(T)sin -^(t-z) ^ X(T)sin -^(t-z) 


(3.33) 


where X(z) is given by (3.23) and 


X(t) = X (t) - (sin 


. 2 TTt .2 7ra.-l/2 


. 4. MClv 

sin 


(3.34) 


Consider, for example, the simple case of an infinite plane containing 
uniformly "loaded" periodic cuts for which fj^(t) = - p^ and £2^^^ “ 

From (3.33) G may be obtained as 


G(z) = ip^ [1-X(z)sin . 


(3.35) 


Using now the relations 


u(x,+0) = •|■[G‘^(x) - G'~(x)] , 


u(x,+0) = - •^[g''‘(x) + G (x)3 , 


(3.36a,b) 


the components of the flux vector on the boundary may be expressed as 
follows; 


u(x,+0) = 


0 , a<|x|<b , 

«,p sin(lTx/2b) 

o 

[sin^ ('rra/2b) -sin^ (xx/ 2b) ] 


, 0<|xl<a , 


u(x,+0) = 


0< X <a 


p sin(Fx/2b) 

— ;° — TT2 - Po ’ s<U|<b (3.37) 

[sin'^(7rx/2b)“sin (7ia/2b)] 


In limit, for b-^, a = finite, (3, 35)-(3. 37) reduce to the following re- 

3 

suits for a plane with a single cut, (-a, a) for which -^ u(x,0) = -p^, 
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(-a<x<a) is the only external disturbance: 


G(z) = ip^[l - 2 (z^-a^) , 


u(x,0) = 


, 2 2.1/2 ’ 

(a -X ) 


0< X <a , 


u(x.4.0) ■ - Po ■ 

(x -a ) 


a< X <00 . 


(3.38) 


In the general case the real and imaginary parts of G(z) give the 
components of the flux vector in the upper half plane and (3.36) that on 
the boundary. g"^(x) and G“(x) may in turn be obtained from (3.33) by using 
(3.24) and the following Plemelj formulas [3]: 


<^(z) = 0^ / 

2771 T . 77 . 

^ S in (t-z) 


, _ Jo, xeLV 

j) (x) (j)“(x) = \ 


(j) (x) + {j)''(x) = 


f (x) , xeL , 


2(j)(x) , xeL' 


1 j f(t)dt 

77i L ■ . 

^ Sin — (t-x) 


xeL , 


(3.39a-c) 


where f(t) is any Holder-continuous function defined on the open interval 
L, L + L’ ^ (-b,b) . 


From (3.13) and (3.33), or more specifically, from (3.37) it is seen 
that the components of the flux vector 9u/3x and 3u/9y have integrable 
singularities at the points of intersection a^, b|^ Mth a 

power of - 1/2 (see also (1,7) for 0^=77). A close examination of (3.13) 
and (3.33) would indicate that 




G(z) = F^(z) + X(z)F 2 (z) 


(3.40) 


where F^ and F2 are holomorphic in the entire plane. By letting z-a = 

X0 

re , for small values of r G(z) may be expressed as 

9 . . e 

.V cos 7T - 1 sxn ir 

n/ \ oa .t3u 1 2 

9x ^ 3y ^ .7T . TTa.1/2 ^2^^^ +0(1) 

/F t^sxn-^] 


O<0<7T 


(3.41) 


For example, in the case of uniformly "loaded" cuts, from (3.35) it follows 

(1) 


nr \ • 5a 


fh . ve.^l/2 f . 0 , . 0 \ 

- (— tan ^) - — (sxn y + 1 cos y) + 0(1) 

0<9<tt-0 (3.42) 


0 

Note that in this simple example u(x,0) = O on L, -5— u(x,0) = - p on L„. 

i dy o 2 

Hence, the boundary conditions around the flux singularity x = a are iden- 
tical to (1.6). It would then be expected that the resulting asjrmptotic 
solutions be the same which may indeed be seen from (1.7) and (3.42). . 


3.3 An Elasticity Problem for a NOnhoraogeneous Plane 

As another example for the direct application of complex potentials 

to mixed boundary value problems consider the following elasticity problem: 

( 2 ) 

Let two linearly elastic isotropic half planes with material constants 
^i**^i ~ ^ for y>0, -°°<x<“, i.e., S and i = 2 for y<0, -“<x<“, i.e,, S ) 

be bonded along the (nonintersecting) line segments = (aj^,bj^) , (k = 1, . .N) 

on the real axis. Let L = K L, be finite and the x and y- components of 


^ ^In (3.41) and (3.42) the terms 0(1) stands for bounded terms and come 
from F.(z) in (3.40). . Note that the second term in (3.40) is of the 
asymptotic form O/r""^'^) + 0(r^/^). 

(2\ 

4 Where ]i. is the shear modulus, k. ~ 3-4v^ for plane strain, and = 
(3-\)i)/(I+Vi) for plane stress, V. being the Poisson’s ratio. 
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the resulting force acting on the half planes at infinity be Q and P, 
respectively. Let L’ be the complement of L on y=0) and -pj^(x) 

and p2(x) the normal tractions acting on the half planes along L’. 


The problem will have to be solved under the following conditions: 


.(» 


, - 
lyy 

. + 

- X cr. 
Ixy 

1 "^(t) - 

lyy 

■ ^ ^Ixy 


[u^(t) + iv^(t)] - [u^Ct) + iv^Ct)] = £^(t) + if2(t) = £(t), tel. 




(3.43a-d) 


where the superscripts + and - refer to limits as yH-0 and y^-0, u^^, v^, 
are the x,y- components of the displacement vector, and 
(j> i) = (x,y) is the stress. 

The simplest method to solve this problem would be the use of complex 
potentials known as Kolosov-Muskhelishvili or Goursat functions. In 
terms of these potentials the stresses and displacements may be expressed 
as follows [2,8-10]: 


\xx ■*' ^kyy " + \(z)] , 


\yy - ^kxx ^^^kxy = ’ 


2y^(Uk + ivj == " z$T,(z) - , 


"k'^k 


\(z) = <j>i;.(z) , T^(z) = , (k-1,2) 


(3.44a-e) 


^^^The input function f(t) may be non-zero Itij for example, thermal stress 
and residual stress problems [15]. 
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*4* — 

where zeS for k =1 and zeS for k = 2. Noting that and are de- 
fined in and $2 '^2 defined in S only, by extending the def- 

inition of into S and into in such a way that they are holomor- 
phic on the unloaded parts of the boundary (i.e., the real axis), one 
could make the following substitution [2]: 

\(z) = - §y.(z) - i^(z) - zf^(z) (3.45) 

■f" — 

where zeS for k = 1 and zeS for k =2. From (3.44) and (3.45) it then 
follows that 


\yy " ^\xy = , (k = 1,2) (3.46) 

^k^^^ - (z-z)l^(z) . (k =1,2) (3.47) 

Substituting from (3.46) into (3.43a) we find 

^l(t) + ^ 2 ^t) = §”(t) + f"(t) , teL + L* , (3.48) 


meaning that §^^(z) + $ 2 ^^) holomorphic in the entire plane Including 
the real axis. Noting that the stress state at infinity vanishes and 
assuming that the rotation at infinity is zero, following [2] , for large 
values of |z| it may he shown that 

+o(l/z) (3.49a,b) 

where, in the usual notation o(l/z)<c/z, c being a positive quantity which 
depends only on jz| and tends to zero as |z|-x». Since §^+$2 holomorphic 
from (3.49) it is clear that 

§j_(z) + = 0 (3.50) 
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In the entire plane. 


Substituting now from (3.46), (3.47), and (3.50) into the boundary 
conditions (3.43b) and (3.43c) (after differentiation) , we Obtain the 
following Riemann-Hilbert problem for the sectionally holomorphic func- 
tion 


r+ 


§ 2 (t) + to § 2 ^^^ "" ’ 

$2^t) - §2^^^ "" p(t) V 


teL 


(3.51) 


where 




(3.52) 


Referring to [3] , the fundamental solution of (3.51) satisfying the related 
homogeneous equations 




X (t) + tox (t) = 0 , teL , 




X (t) - X (t) = 0 , teL’ 

k ’ \ 


N 

L = E L 
1 


(3.53) 


may be exp re s s ed as 


N 


X(z) = (z-bj^)^(z-a,^)Ac 




\ 27ri \ 2 ^ 2iT 






(3.54) 


V7here and (k = 1,. .,N) are again arbitrary (positive, zero, or neg- 


ative) integers In the present problem the singular points a^ and 
b^ correspond to the ends points of interface cracks. Consequently, at 
these points the stresses and displacement derivatives will have an in- 
tegrable singularity giving = 0 
then consider the particular branch of X(z) which is single valued in 
the plane cut along L and for which 

lim z^X(z) = 1 . (3.55) 

\z\-^ 


, = -1 , (k ■- 1 , . , ,N) . Here , we will 




Dividing by (3.51) by X (t) and using (3.53) it is found that 
4j(t) + - _ h(t) 

\ v/x-N/ 


(- 


X(t) 


•) 


X(t)- 

«>o(t) 


x'^(t) 


^x(t)> - ^ X(t)^ - 


teL 


teL’ 


(3.56) 


Noting that the stress state vanishes at infinity, the general solution 
of (3.56) vanishing at infinity becomes 


N-1 


, 

X(z) 27TX ^ (t_z)x^(t) L' (t-z)x'^(t) 0 


^ / _Mt)d^ + 77^ / — -f-s c^z” , (3.57) 


where c^,..,c^ are arbitrary constants. From (3.49b), (3.55), and (3.57) 
it may be shown that 


“N-1 


CH-iP 


2tt 


(3.58) 


The condition of single- valuedness of displacements gives the remain- 
ing N-1 constants, c^,. .,Cj^_ 2 . We recall that in deriving (3,51) the 


^ ^In most physical problems and are such that -l<Re(a]^, 3|^) <1 . Even 
though in literature one finds this as a mathematical condition, clearly 
Ajj_ and Bjhmust be determined from the physics of the problem. 
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continuity condition (3.43c) was used in differentiated form. Thus, this 
condition requires that 


- 9 + - 

\ 3^ [(ui+ivj) - (u^+iVj) ]dx = - f(b^) 


(k = 1,..,N-X) 


(3.59) 


or 


a, 1 1+K- I+K 2 , 1 *^2 


(k = 1, . . ,N-1) . 


(3.60) 


In (3.60) and elsewhere, to obtain the boundary values of the Cauchy 
integrals the following general Plemelj formulas may be used 

FCz) . ^ 

2'rr-x ^ t-z ’ 

u 


f‘^(x) - F“(x) 


F (x) + F“(x) 


g(x) 

.0 , 


xeL 

xeL’ 


Fx ^ t“X ; 

2F(x) , xeL’ 


xeL 


(3.61a-c) 


where L’ is the complement of L (on an infinite line or on any closed con- 
tour in the complex plane) . Thus, (3,57) with (3.58), (3,60), (3,50) ^ 
(3.45) and (3.44) gives the complete solution of the problem. For example 
if If = (-a,a) (bonding along a single segment) , p(t) = 0, and f(t) = 0, 
the solution becomes 

i(logW)/2F 


$i (z) = - fo(2) 


Q+IP 




2i , 2 2,1/2 ^z-a 

(z -a ) 


(3,62) 


where OJ is given by (3.52). 
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Observing that the general solution of the problem is of the form 


= F^(z) + F2 (z)X(z) , (3.63) 

as in the previous example, the asymptotic behavior of the stresses and 
displacements around the points of singularity may easily be investigated 
(see [16] for details). If the positive constant 0 ) appearing in (3.51) 
and defined by (3.52) is not unity, from the behavior of the fundamental 
function X(z) given by (3.54) it is clear that the stresses and displace- 
ment derivatives will have a typical oscillating singularity around the 
end points of the branch cuts which is of the following form: 

aj_j(r,6) = 0)cos(Y log j) + g_(0)sin(Y log j)] + 0(1) , 

Y (log w)/ 27T , (i,j = x,y) , (v«t) (3.64) 

where r,0 are the polar coordinates around the singular point, f^^ and 
g.j are bounded functions, and the term 0(1) again comes from the analy- 
tic function fj^(z) in (3.63). 

4. REDUCTION TO SINGULAR INTEGRAL EQUATIONS 

As pointed out in the previous section, in order to obtain the 
correct behavior of the singularities and also in most cases to find 
a simple closed form solution of a given mixed boundary value problem, 
whenever possible it is always preferable to use the complex potentials 
and the related complex function theory. However, the technique has its 
limitations. First, the particular problem may not admit complex poten- 
tials. Secondly, the successful application of the technique is severely 
limited to certain geometries. Finally, there are certain types of 
boundary conditions which: would make the direct use of complex potentials 
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xrawffy,imagg»iigff3gw*^:r^^ 


extremely difficult if not impossible (e.g., the boundary conditions 
described by differential operators containing the unknown function as 
well as its derivatives . In such cumbersome cases, the method which 
is perhaps the most general and the easiest to apply is the reduction 
of the problem to singular integral equations either by using a Green's 
function formulation or by formulating the problem first in terms of 
multiple series or multiple integral equations. Reduction of the bound- 
ary conditions to integral equation is always possible. The main prob- 
lems here are the selection of the appropriate auxiliary function (i.e., 
the new unknown function defined on the boundary) and the proper separa- 
tion of the dominant parts of the kernel for the correct study of the 
singular behavior of the solution. In this section this important step 
of reducing the mixed boundary value problem to a singular integral 
equation will be discussed by considering some typical examples, and 
some general remarkswill be made regarding the nature of the kernel , 
and the solution. 

A.l Reduction of Dual Series Equations to 
Singular Integral Equations 

Consider the typical simple mixed boundary value problem described 
by (2. 1-2. 3) and formulated by the dual series equations (2.6). At the 
generality that the problem is stated, it is not very fruitful to pursue 
a technique based on the operational methods to solve the problem and, 
as stated before, because of the importance of existing singularities, 
a direct numerical solution is nearly useless. To reduce the problem to 
an integral equation, the first step is the definition or selection of an 

^^^Eor this type of mixed boundary conditions even the Simplest problems 
such as that in potential theory for a half plane do not seem to have 
closed form solutions (see the following section) . 
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appropriate auxiliary function. In this problem let this function be 

(f)(e) = -^ u(l,0) , 0eL^ + L2 = (0,2tt) . (4.1) 

Note that 

4)(e) = 

From (2.5) and (4.1) it follows that 
nA = - ~ /j ,T <})(t)sinnt dt , 

T1 7T * ^2 

nB^ = ^ A (J)(t)cosnt dt . (4.3) 

n TT 1']^+L2 

Equations (4.3) and (2.5) give the solution once the function (j)(t) and 

the constant A are determined, 
o 

In applying the technique described in this section, in order to 
circumvent the difficulties arising from the divergent series or integrals 
giving the kernels, for analytical convenience the boundary condition form- 
ing the basis of the integral equation will always be expressed in limit 
form^^\ Thus, the boundary condition (2.6a) may be expressed as 

lim E n r” ^(A cosn9+B sln-nG) - g(0) , 0£L, (4.4) 

r^l-0 1 ^ 

Substituting now from (4.3) into (4.4), observing that for r<l the related 

Infinite series will be uniformly convergent, and hence changing the order 

^^^The exception here is, of course, the case in which the related series 
or integrals are uniformly convergent giving bounded kernels . In that case 
the limit may be put under summation or integration sign before evaluating 
the kernels. 



0eL 


2 

0EL 


1 


(4.2) 






of integration and summation, it is found that 

lim A 4>(t)dt E r” ^sinn (t-^0) = - g(9) , 6eL, . (4.5) 

T^l-0 ^ 1 ^ 

First performing the sum as 


“ n-1 _, 1 “ r. iZxU , -iZxUi sin z 

E r sinnz = -r^ E Rpe ) - (re ) J = z — 

1 1 l+r'^-2r cos z 


(4.6) 


and then going to limit, (4.5) becomes 

^ /j.^4>(t)cot dt = - g(9) - •'■L^f'(t)cot dt, 9eL^ (4.7) 

giving an integral equation to determine (J)(6). 


It should be observed that (4.7) is a typical singular integral equa- 


tion of the following form 


? \ rj,^k(0,t)*(t)dt = g^ce) , eeLj 


(4.8) 


where the kernel 


1 ^ .-N 1 t-6 11 

k(x,t) . ^ cot — - -jre 


(4.9) 


is bounded everywhere on Lj^ (including the end points a j^,b^, i = 1,,. ,N), 
and the known function g^ is the right hand side of (4.7). Since the 
fundamental solution (or the behavior of the singularity) depends only on 
the dominant part of the Integral equation, for the purpose of obtaining 
this solution tentatively expressing (4.8) as 




(4.10) 


defining a sectlonally holomorphlc function by 


“ 4r A ^ ' 
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and using the Flemelj formulas 


f‘*‘( 0) - F"(0) = (J)(0) , 0eL^ , 

F-^(e) + F-(6) = i ^ dt , 
one obtains the following Riemann-Hilbert problem: 


(4.12a,b) 


F (0) - F (0) = 0 


0eL. 


F'^(0) + F"(0) = - igj(0) , 


0eL, 


(4.13) 


N 


Again, following [3] and observing that L^, = E funda- 

mental solution of (4.13) satisfying the related homogeneous Riemann- 
Hilbert problem is found to be (see (3. 5)- (3. 10)) 


N 


X(z) = n(z-a.) )l/2+efe 


(4.14) 


1‘ "k' ‘'k" 

where the arbitrary integers and 3^^ (k = 1,..,N) have to be selected in 
such a way that the solution is consistent with the expected physical be- 
havior. In this problem since at a^ and bj^ the "flux" has an integrable 
singularity, Oj^ = 0 and 3j^ = -1, k = 1,..,N. Referring, for example, to 
the solution of (3.5) as given by (3.12) and noting that X^(0) + X (0) = 0, 
OeLj^, from (4.12a) it is seen that ({>(0)'“X"^(0) , 0eLj^. We now define the 
fundamental function, w(0) of the singular integral equation (4.10) by 


N 

w(0) - qi [(e-a;^)(b^-e)]"^'^^ = (-1)^'V(0) 


(4.15) 


Thus, the solution of (4.10) or (4.7) will be of the following form 


(j)(0) = w(0)p(0) 


(4.16) 


where p is an unknown bounded function. 


b 






Referring again to (3.12) (and [3]), it is seen that the solution 
of (4.13) vanishing at infinity (and hence <{)(0) obtained from (4.12a)) 
will contain N arbitrary constants. These constants are determined from 
the single-valuedness conditions of u(rj0). We recall that in deriving 
(4.7) the boundary condition u(l,0) = f(0), 0eL2i was used in differ- 
entiated form. Therefore (J)(0) must satisfy the following conditions: 
b. 


(J)(0)d0 = f(b^) - f(a^) , (k » 1,..,N) . 

k 


(4.17) 


To complete the solution of the problem the constant in (2.5) must 
be determined. From (2.5), (2.3), and (4.1) it is seen that 

1 1 N ^k+1 

A = -h A XT ua,e)d0> ? V f<0)d0 

o 2Tr Lj^+L 2 2ir 1 bj^ 


+ ^ d0[f(a^) + (|)(t)dt] , = ap . 

3 3 


(4.18) 


Using the Dlrichlet transformation 


dy F(x,y,s)f(s)ds = f(s)ds 7* F(x,y,s)dy 
a a a » 


(4.19) 


equation (4 . 18) becomes 

1 W 


K " 2U i 


f (0)d0 - /7 t(})(t)dt + bj^f(b^) - ay.t{a^)] , 

■ k 


S+i " • 


(4.20) 


The solution of (4. 13) is given by (3. 12) with gj^ - 2fj^, f^ = 0, and 
N - 1 as the degree of the polynomial P^. Thus, using the Plemelj formulas 
(4.12), the solution of (4.7) or (4.8) maybe expressed as 
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( 4 . 21 ) 


N-1 8 (t)dt 

(j)(0) = 2w(6) E c e"" - /- 

^ O n IT Lj^ (t-0)w(t) 


+ /. 


dt 


IT 'Lj (t-e)w(t) • 


Considering now the definition (4.16), it is found that 


N-1 


p(9) - “(®»s)p(s)ds = 2 E c^0^ - ^ / 


gQ(t)dt 


o n 


TT Lj^ (t-0)w(t) ’ 


0eLj (4.22) 

where 


n(0,s) 



k(t,s)v7(s) , 
(t-0)V7(t) 


(4.23) 


Equation (4.22) may be treated as a Fredholm integral equation [3,5] giving 
the solution of the problem with (4.17) and (4.20). 


It should be remarked that any singular integral equation of the form 
(4.8) may be "regularized" and reduced to a Fredhom type integral equation 
by following the foregoing procedure. However, for the solution of (4.8) 
a considerably simpler numerical technique will be described later in this 
article. It should also be noted that the particular mixed boundary value 
problem considered in this section can be solved in closed form. First 
consider the case of geometric symmetry with respect to the line 0 = 0. 
That is, let us assume that 


m 


N = 2m 


. h - ? ’ ”lk ' <vV • 


m 


“i " ^ '•ik ’ “ik -/-'>k--"k> • 


(4,24) 


Then the problem may be reduced to a singular integral equation with dom- 
inant part (i.e., Cauchy kernel) only. For this the problem is first 


decomposed into symmetric (i.e., u(r,6) = u(r,-0)), and antisymmetric, 
(i.e., u(r,6) = - u(r,-6)) parts by separating the input functions into 
even and odd components. In the symmetric problem g(0) = g(-0), f(0) = 
f(-0), (f)(0) = - cf)(-0), and (4.7) may be expressed as 


1 H -^a^ 4>(t)[cot + cot -^]dt = g (0) = - g(0) 

i. k 

- /jj f *(t) [cot + cot -5~]dt , 0eMj^ 


(4.25) 


where 1^2 is the complement of on O^O^TT. Or, from 
_ sin 2 b 

cot(a-b)+cot(a+b) = 2{ . _ }/(cos2b- cos 2a) 

SXH Z3- 


(4.26) 


it is found that 


1 r cj)(t) sin t 
■ J ' ■ Q ' ■ dll 

7T COS 6-COS t 


go(e) 


0d^n 


(4.27) 


If we now define 


cos t = a , cj)(t) = t|J(a) , cos 0 = 3 , 8^(6) = G(3) , 

m 

teM^^ aer , r = 2 Tk (4.28) 

equation (4.27) becomes 

i f 

Similarly for the antisymmetric problem g(6) - - g(-6) , f(0) = - f(-6) 
(f)(0) = ^(-0) and (4.7) becomes 

In ^ ^]dt = g^(0) = - g(0) 

- ^ f ' (t) [cot - cot -^—Idt , 0eMj^ , (4.30) 
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with the substitutions (4.28) giving 
j. f i|/(g)\ da / G(g) 


\/l-a 




3er 


(4.31) 


In both cases the solution is given by (4.21) with k(t,s) = 0 and approp- 
riate changes in notation, m integration constants arising from the solu- 
tion are again determined from the single-valuedness conditions (4.17). 

In the symmetric case is given by (4.20) (with, again the appropriate 
changes) and in the antis^Tnmetric case it is zero. 


For example, if 

m - 1, = ( 0 j^, 02 ) , a^ = cos 0j , a 2 = cos 02 , 

O<0^<02<tt , g( 0 ) = g(-e) , f(0) = 0 , (4.32) 


it may be shown that 


j ^ a3^G(3)(3-a2)^^^(aj^-3)^^^ 

[(a-a2)(a^-a)]^/^ 


ijj(a) = 


0 , (-l<a<a2 , aj^<a<l) 


(a2<a<a^^) , 


(4.33) 


where c and A are determined from 
o o 


a. 


j. 1 i|j(a)da ^ Q 

A-a^ 


dr 


1 

A = - — /q t(j)(t)dt 
o ^ 


(4.34) 


Furthermore, if g(0) — a^ = constant, ip(a) becomes 


- 1/2 - 1/2 ^ l "*^2 

ip(a) = ^ (a^-a) - a + b^) , 


l+a^ 


b = 
o 


a„~a, 2+a,+a„ 

n (4i^. w - 


l+a. 


2 
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2(aj-oij) 







I' 




n (c?,k) /q 2v. 2v1/2, 2 2.1/2 

(1+ax )(l-x ) (1-k X ) 

where K(k) is the complete elliptic integral of the first kind. 


(4.35) 


Consider now the general case (4.7),i.e. 


W dt = 8 o(Q) » 9eLj 


(4.36) 


where 8^(9) is the right-hand side of (4.7) and = Z has no symmetry. 
Let the ’origin is selected in such a way that 0 = 7T is not on Lj . Defining 


t 0 

tan j = s , tan j = p , 4>(t) = i]^(s) , 8^(6) = G(p) , 


teL, sEr, , r, = z r., 

1 1 ’ 1 1 Ik 


(4.37) 


equation (4.36) may easily be expressed as 


i / ^.(s) ^ - G(p) , 


(4.38) 


1 / iki ^ ^ g. 

TT p S-p 


(4.39) 


K = i / 


1 r s^Ks) 


(4.40) 


The singular integral equation (4. 39) may be solved in a straightforward 
manner with (4.40) accounting for the additional constant K. Again, the 
solution of (4.39) will contain N arbitrary constants which may be de- 
termined from (4.17) and (4.20) gives the constant A^. 


4.2 An Example on Triple Series Equations 

Consider now the somewhat more general mixed boundary value problem 
for the unit circle defined by (2.10) and formulated by the triple series 
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equations (2.11). In this problem let the auxiliary function be 


( 1 ) 

(})(0) = u(l,0) , 0 eLj, + L 2 + L 3 = L = (0,2it) . (4.A1) 

From (2.5) it then follows that 

nA = — /_ (f)(t)cosnt dt , 
nTTL^ ’ 

nB^ = ^ (j)(t)slnnt dt (4.42) 

where 4>(t) = fj^(t) is known on Lj^. After differentiating (4.11b) and, 
again for analytical expediency expressing them in limiting form, the 
remaining boundary conditions (4.11b and c) may be written as 

6 eL 2 

.OeLg (4.43) 

Substituting now from (4.42) into (4.43), using (4.6) and [17] 

log (2 1 sin j\) , (4.44) 

we obtaj.n 

(|)(t)cot dt = £2^0) , 0St2 ’ 

h. ' 

— — /j^4>(t)log(2 I sin “Y"! )(it + h 2 <j>( 0 ) = ~ h^A^ jGeL^. (4.45a,b) 

^^^See the general remarks and the broad guidelines at the end of Section 4 
regarding the selection of the auxiliary (i.e., the new unknown) function. 


00 

H — cos nz = - 

1 ^ 


“ n 

lim S n r (-A sinn0+ B cos n 6 ) = fl(0) , 

r^l -0 1 " ^ 

CO 

lim {h,A + h. 2 r (A cos n0 + B sinn9) 
r->l-0 1 “ ll ” n 

+ h^E nr” (A cos n0 + B sin n0) = f „ ( 6 ) 
2i n n ' 3 
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Just a superficial observation would indicate that the two integral 
equations in (4.45) are of entirely different character. Hence, near 
and at the end points of L 2 and one would expect the function (j)(t) 
behave quite differently. For closer examination let us assume that 


'fi(t) = 4>2(t) , teL 2 , = ^ 2 ^*^) » *^^^3 


(4.46) 


Equations (4.45) may then be expressed as 




7T 


= £ 2 ^ 0 ) - /j_^f^(t)cot dt = gj^(Q) , 0eL2 , 

Vi 

/l <f)j^(t)log(2|sin -^“hdt - “ (f)2 (t) log(2 | sin -“| ) 


n. 0 

+ h2<f)2(0) = ^ 3 ( 6 ) “ IT ^i(t)log( 2 |sin -^-1 )dt = 62 ( 0 ) j 


0£L, 


(4.47a,b) 


Separating the dominant parts of the kernels it is seen that (4.47) is of 
the following general form: 

1 (t) 

¥ kj2(9,t)<|.2(t) dt , 


0eL 


2 ’ 


n, 

h2(f)2(0) “ ^ = §2^0) “■ 


- fj^ k22(9.t)(j)2(t)dt , ©eL^ 


(4.48a,b) 


where the kernels k^^( 6 ,t), (i,j = 1 , 2 ) are bounded in their respective 
domains . 
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It Is clear that (4.48a) is a typical singular integral equation 
of the general form (4.8), it has a fundamental function of the form 

(4.15), and if L 2 - ^ ^2k* ^2k ~ ^^k’^k^ ’ ^k^^^k'^^k+l ’ solution will 
contain N arbitrary constants. On the other hand, (4.48b) is an integral 
equation of the second kind with a weakly singular but square integrable 
kernel. Hence, its solution is bounded everywhere on L^, Including the 
end points (as they are approached from L^) and is uniquely determined 
without any reference to any additional conditions. The coupling of the 
tv;o integral equations is through Fredholm kernels. Therefore the basic 
singular behavior of the solution will be unaf fected by coupling. 

From (4.10) it is clear that if hj^ or h 2 is zero the problem reduces to 
that considered in the previous section. This can also be seen from (4.47) 
For hj^ = 0 this is quite clear. For h, = 0 differentiating (4.47b) and 
observing that 

- ^ log(2lsin -^^1) = j cot (4.49) 

(4.47) is seen to reduce to a simple singular integral equation defined on 
^2 ^ ^3* worthwhile to reemphasize that an integral equation of the 

second kind with a logarithmic kernel is basically a Fredholm integral 
equation and has a bounded solution. However, if the integral equation 
is of the first kind and has a logarithmic kernel, then it is equivalent to 
a singular Integral equation with a simple Cauchy type kernel. Also, the 
dominant equation 

A(b(x) + B /j^{f)(t)log j t“x| dt = f(x) , xeL (4.50) 

As will be pointed out later in this article, if the coupling is through 
generalized Cauchy kernels, then the singular behavior of the solution will 
be affected . 
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has apparently no closed form solution. Even though at the end points 
of L the solution of (4.50) is bounded, it may easily be shown that its 
derivative has a logarithmic singularity. This can be seen by differ- 
entiating and modifying (4.50) as follows: 




f f KtMCsI <jt + f 4i(x) <• 

A L t-x A 


■^L t-x 


dt ^ f * (x) 


xeL 


(4.51) 


which, in the neighborhood of a typical end point x<c becomes 


c|j'(x) = j (f'(c)log|c-xl + G(x) , (x^c) 


(4.52) 


where G(c) is bounded. 

N 

If L- = , to complete the solution of the problem the arbitrary 

^ ^ K. 

constants c ,c,,...c„ , arising from the solution of the singular Integral 
o 1 N-1 

equation (4.47a) and A^ must be determined. Observing that the boundary 
condition u(l, 0 ) = f 2 ( 6 ), satisfied only in differentiated 

form and 


u(l,0) = A - / (|)(t)log(2lsin -^j)dt , 

■ L 


(4,53) 


4 >(t) must satisfy the following single-valuedness conditions: 


t-e. 


A^ " ~ / ^(t)log( 2 | sin ■ ^ - - |)dt - ^ 


O 7T 


, k = 1,. .,N 


(4.54) 


where 0,^ is any convenient point on L 2 j^. Tbius, equations (4.54) with the 


flux equilibrium condition 
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fj(t)dt + (j)2(t)dt = 0 (4.55) 

provide N + 1 algebraic equations to determine and A^. 


4.3 Reduction of Multiple Integral Equations 

To demonstrate the technique of reducing a system of multiple in- 
tegral equations to that of singular integral equations consider the 
somewhat general mixed boundary value problem in potential theory which 
is defined by (2.21) for the half plane y>0. The solution is expressed 
by (2.17) in terms of the unkno^ra function A(a) which is to be determined 
from the multiple integral equations (2.22). In this problem let the 
normal flux 


u(x,o) = (|)(x) = -|a[A(a)e ^°^^da , -<»<x<°° (4.56) 

be selected as the new unknown function. If 4>(^) is known, the Fourier 
inversion 

I !*/ N 1 (4.57) 

-la|A(a) = -^ /^ (j)(t)e dt 

with (2.17) would give the complete solution. From (2.21) it is seen 
that c|)(x) = ^l(^) i® kno^TO on and is unknown on L 2 and L^. Again for 
reasons of uniform convergence writing the boundary conditions (2.21b) 
and (2,21c) in limit form, differentiating (2.21b) and using (4.57) we find 


T - fi ~“(y+ix)j A/-. N iottj.. 

lim [-^ / e ^da (|)(t)e dt 

jrH-0 ° 


i rO a(y-ix)j^ 

- 2 ^ da (j)(t)e dt 


f^(x) , xeL 2 


1 , 1 • r "Ct(y+ix) da r°° 

^ h^ lim [- f^e ~ c|)(t)e dt 


y^O 


41 


(4 .58a,b) 


+ !l 4 r C «t)e‘«dtj + hj^Cx) = fjCx) 


a -oo 


xeL 


3 • 


Changing the order of integrations and evaluating the inner integrals, 
(4.58a) becomes 


T 00 1 1 

■ 7 ^ lim / <p(t) [ — TTz — •. ' , ~ /Z — \ 

,q y-x(t-x) y+i(t-x) 


= t'(x) . 

IT -00 t-x 2 ’ 


xeL, 


(4.59) 


Similarly, the first term on the left hand side of (4.58b) may be ex- 
pressed as 

h. -cxy 

lim / (|)(t)dt/ — cosa(t-x)da 

tt -00 / o a 


7T 


-ay 

00 ^ 00 A 

lim (j>(t)dt{/_^ + log y 

3TH-0 " 


o a 


log[y^ + (t-x)^]^'^^} = (f)(t)loglt-x|dt 


1 


(4.60) 


where the following condition of flux equilibrium is used to eliminate 


the divergent terras: 


. (1) 


/ (|)(t)dt = 0 


(4.61) 


Now observing that - ("*°°j°°) > 4’(b) = f j(t) , tcLj^, and 

defining 


(|)(t) - (j)^(t) , teL2 , 4)(t) = (P^(t) , teL^ 


(4.62) 


Note that if (4.61) is not satisfied and if u is zero at infinity (as 
assumed in the present problem) then u(x,0) will not be bounded and will 
tend to infinity as log y, y->“0. 
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the Integral equations (4.58) may be expressed as 


1 1 
ifj dt+^f. 

TT t-X 


fjCt) 


1 

TT t-X “ “ ^ 2 ^^^ “ ? t-x 


dt 


xEL^ 


h h 

+ ~ (j)2(t)logit-xidt + -Iff (|)j^(t)loglt-xldt 

= fi(t)logl t-x|dt , xeL^ . (4.63a,b) 

In structure the system of equations (4.63) is identical to (4.47) or (4.48) 

N 

If L„ is finite with L„ = EL_ , , . = (a, ,b,), the solution of (4.63a) will 


2j’ ^2j 


j’ 3' 


again contain N arbitrary constants which may be obtained from (4.61) and 
the following N-1 single-valuedness conditions: 


i /” (|)(t)log 


t-a 


j±l 


t-bj 


dt = f2<»J+l' - 


j - 


(4.64) 


In the type of problems under consideration the case of non-zero re- 
sultant flux may be particularly important. From (4.60) it is clear that 
if (4.61) is not satisfied the derivation leading to (4.63) is not valid. 


Let us now assume that is finite 


( 1 ) 


and 


/ (j)(t)dt = Q 


(4.65) 


^since u(x,0) = f 2 (x), xeL 2 is assumed to be bounded) u(x,y)~log y as yx». 
lowever, the components of the flux vector may still be expressed in terms 
)f an unknown function A (a) as follows: 

^ u(x,y) = - / iaA(a)e ^*^^da , 


Sx 


^^^If L 2 is infinite, by a proper superposition the problem may always be 
reduced’to that in which (4.611 is satisfied. 
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(4.66a,b) 


u(x,y) = - l2 |ot|A(a)e~^^^l“^“’^da . 

Again, defining the auxiliary function <{)(x) by (4.56) and substituting 
from (4.57) into (4.66a) we obtain 


3 i T 4 r -Oi(y+ix) j r°° 

u(x,0) = — lim 1/ e 'da / <P(t)e dt 

(7X ZTT _vii^ O — ^ 

3TH-0 


.0 a(y-ix) ... 

/ e da / (p(t)e dt] 


i . r°° j. /_N 2i(t-x) 

= ^ L ^(t) -0 T dt 


2 

' 3TH-0 


9 9 

y +(t-x)“ 


*(t) ^ 

TT t-x 


(_oo<x<oo) 


(4.67) 


from vjhich it follows that 


u(x,0) = ^2^^o^ + ^^t)log|-^ 3 ;p|dt , (-oo<x<oo) 


(4.68) 


where x^ is any point on L 2 . Thus, the integral equation (4.63a) is still 
valid and, with the definitions (4.62) and using (4.68), (4.63b) will have 
to be replaced by 

h^+jCK) + ^ j- <t,^(tnog ^ /j_j.^(t)iog -g- dt 


fj(x) hjfjCx^) ^ dt , xeLj 

1 O 


(4.69) 


the N integration constants arising from the solution of (4.63a) are determ- 
ined from (4.64) (which follows now from (4.68) and from the fact that 
u(x,0) = f 2 (x), xeL 2 ) and (4.65). Note that under the stated conditions 
of the problem (4.63b) and (4.69) are always nonhomogeneous. Also note 


that if I 2 - problem is reduced to 


1i2(j)2(x) + — /t. 4>9(t)logjt-x]dt = £(x) , xeL, 


IT 


(4.70) 


v;here f(x) is a known function. !5quation (4.70) has no closed form solu- 
tion. In this case it should be noted that u(x^,0) = f2(x^) which appear 
in (4.69) is an unknown constant and is determined by the flux equilib- 
rium condition (4.65) where x^ now is an arbitrary point on the real axis. 
On the other hand if - 0 and L2 7^ 0, the problem is formulated by 
(4.63a) which may easily be solved in closed form. 


As pointed out previously the solution of the system of singular 
integral equations (4.63) will be of the following form 

N 


cj)j(x) = w(x)Fj^(x) , w(x) = f~| [ (x-a^) (b^^-x) ] 
4>2(x) = 


- 1/2 


xeL^ 


(4.71) 


where and F2 are bounded functions and at the end points of 4>2 
bounded and has a behavior similar to (4.52) whereas at the end points 
of L2 is singular. 

In order to have a better understanding of the flux distribution 
^(x) on the boundary it tnay be worthwhile to relate the foregoing results 
to a simple physical problem. Consider, for example, the anti-plane 
shear problem for a s3nnmetrically loaded infinite medium shown in Figure 2. 
Here u(x,y) is the z- component of the displacement vector* <j)(x) = CJy^(x,0)/u 
represents the traction at y = 0, and pQ is the resultant force 5.n z- direc- 
tion, p being the shear modulus of the medium. The part of the x-^axis 
L = El,, corresponds to a series of cracks on which the surface traction 
(j)(x) = f^(x) = cTy^/p ' is specified. L2 is clearly the uncut portion of 
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the real axis on which (in this problem) u(x,0) = ^2^^^ ~ Pl^ysically 
the part on which 

hj^u(x,0) + h, u(x,0) = (4.72) 

corresponds to a series of cracks the surfaces of v/hich are joined through 
an elastic adhesive layer. If the thickness of the adhesive is sufficiently 
small it may be represented by a one-dimensional spring and may be modeled 
by (4.72) with 

hi - y^/h , h^ = - y , ^3^^) = 0 . (4.73) 

where y and h are, respectively, the shear modulus and the thickness of 

3 . 

the adhesive. Thus, the spring force P^u/h will act as a traction on 
the crack surface along L^. Physically then the plane is cut along 
Li + L^, which means that the stress a^^(x,0) = yi|)(x) must have the ex- 
pected square root singularity at the end points of ^2 (as x approaches 
these points from I- 2 ) • Note that whether the end point belongs to the 
intersection of L 2 and Li or L 2 and L^, this singular behavior will remain 
unchanged. On the other hand in the adhesive layer (i.e. , on L^) the 
shear stress will be bounded everywhere. Including the end points. Further 
applications of the technique to the mechanics of bonded joints may be 
found in [18] and [19]. 

4.4 Reduction of Multiple Series-Multiple Integral Equations 

In some mechanics problems because of the geometry of the domain 
the separation of variables technique may lead to a formulation in which 
some of the unknown functions are expressed as series of eigenfunctions 
having a set of undetermined coefficients and some as inversion integrals 

46. . : 


involving certain undetermined functions. Substitution of these ex- 
pressions into the mixed boundary conditions invariably gives rise 
to a system of equations involving both multiple series and multiple 
integral equations. Once the auxiliary functions are properly selec- 
ted, reduction of this type of problems to a system of singular integ- 
ral equations is also rather straightforward. In this section we again 
select a highly representative and a relatively simple example to dem- 
onstrate the technique. 


Consider the "load transfer" problem shown in Figure 3. The figure 
describes either a coupling in which the torque is transmitted from the 
shaft 1 to the shaft 3 through the sleeve 2, or a gear or pulley in which 
the external torque acting on the sleeve 2 is not zero. In this problem 
it is assumed that the shafts are made of the same material with shear 
modulus and the sleeve has the shear modulus Pp* further assumed 

that (in addition to axial symmetry) x = 0 is a plane of symmetry (or 
antisymmetry) with respect to the external loads and the geometry of the 
problem. Thus, in both media the circumferential component of the dis- 
placement u^, (i = 1,2) is the only non-zero displacement which satisfies 
the following differential equation: 


9^u. , 3u. u, 3^u. 

3x2 


(r<a ; i = 1 , 


a<r<R ; i = 2) (4.74) 


the nonvanishing stress components are given by 
3u . 


^ir0 


P.(-^ 

X 3r 


u , 


0.0 1 ' • 
X0X X 


3u . 

X 

3x 


(1 = 1 , 2 ) 


(4.75) 
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For the sake of generality at this point it is assumed that the end clear- 
ance 2b j is not zero and 2c<2b (Figure 3). Because of symmetry it is suffi- 
cient to consider one half (x>0) of the composite medium only. 


In this problem the solution of (4.74) must be found subject to the 
following conditions: 


g 

[u 2 (a + 0,x) 

- Uj^(a - 0,x)] = f (x) 

, bj<x<c , 

(a) 


(a,x) , bj^<x<c , 


(b) 

'llreCa.x) - 0 

5 x>c , 


(c) 

a2r0(a,x) = 0 

, 0<x<bj^ , c<x<b , 


(d) 


, » 


(e) 


, a<r<R , 


(f) 

Q’2r0(R»^V u^XR.x) = 0 , 

, 0<x<b , 

(g) 

I^UjCr.O) . 0 

(symmetric case) 

! 

1 

> a<r<R 

(h) 

U2(r,0) = 0 

(antisymmetric case)^ 

1 


27Ta^p(x)dx = 
b. 

T , 


(i) (4 


where the functions f(x), CT^(x) and the constant T are known. Referring 
to Figure 3, the, following symmetry conditions must be satisfied: 

For the symmetric problem: 

ii. (r,x) = u,(r,-x) , a. p|(r,x) = a. Q(r,-x) , ± - 1^2 , 
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9 


(4.77) 


f(x) = - f(-x) , = T 2 , cf^(x) = cr^(-x) 


For the antisymmetric problem: 


u^, (r,x) = - u^, (r,-x) , a^^Q(r,x) = - a^,^p^(r,-x) , i = 1,2 , 


ir9' 


f(x) = f(-x) , Tj = - T 2 , “ ~ 


(4,78) 


Considering (4,77) and (4.78) the solution of (4.74) satisfying (4.76e,f 
and h) may easily be expressed as 


Uj^(r,x) 


2 r°° 


^ A(a)I.j^(ar)cosa(x-bj)da , x>b^ , 0<r<a 

CO 

u„(r,x) = Z[B K, (a r) + C I, (a r)]/ } , 0<xb , 

2 ’- j nln nln 1 sxna^x ' - 


a<r<R , 


(4.79a,b) 


where from (4.75) and (4.76f) it is seen that 


= 7rn/b for symmetric problem 


a = (2n - l)Tr/2b for antisymmetric problem 
n 


(4.80a,b) 


In (4.79b) the upper and lower terms in the series stand for symmetric 
and antisymmetric problems respectively. Substituting now from (4.79) 
into the mixed boundary conditions (4.76) (a-d) we find 


- sina xl 
n 


lim Z [B^K (a r) + Cl, (a r)] 
ir\ xi riX Ti nx n 

r->-a+0 1 cosa x 

n 


2 


+ lim — / aA(a)I, (ar)sina(x-b.)da = f (x) , b,<x<c , (a) 

r->-a-0 ° ^ 


i cosa XI 

# _ I 

sxna X 
n 
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2y 


1 


TT 


A(a)al2(oia)cosa(x-bj^)da = 0 , bj^<x<c 


(b) 


2y 


1 


IT 


A(a)al 2 (ota)cosa(x-bj)da = 0 , x>c , 


(c) 


i cosct xj 

" > = 0 , 

sinot X 
n 


0<x<bj^ , c<x<b , 


(d) (4.81) 


Note that one set of constants or may be eliminated by using the con- 


dition (4.76g). For example. 


( 1 ) 


B K,(a R) + C I- (a R) = 0 for u^(R,x) = 0 
n i n n i n i. 

/. 


\ 


- B K (a R) + c I (a R) = cr 

n 4 n n ^ n oy 2 ° 


cosoi X 


o(^)|sina\j 


dx , 


for 02^.0 (R,x) = cTq(x) , 


(4.82a,b) 


With (4.82), (4.81) provides a system of dual series-dual integral equations 

to determine the set of unknown constants B (or C ) and the unknown func- 

n n 

tion A(a) . 


In this problem, the contact stress 


aj^g(a,x) = a 2 j.g(a,x) = p(x) 


(4.83) 


suggests itself as being the most appropriate auxiliary function. Thus, 

from (4.82) and the expressions 

2y, 

^ 1 00 

cr ' a(a,x) = / A(a)aI„(aa)cosa(x-b, )da = p(x) 

lr0 TT o 2 1 


^^^In (4.82b) a (x) = 0 being the practical case of coupling. 
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j, . • 1 


(b^<x<°°) 


(J2^g(a.x) = M2E[- B^K^Ca^a) + C^I^Ca^a^a^ 


cosa X 
n 

sinct X 
n 


= p(x) 


(0<x<b) , 


(4.84a,b) 


evaluating A(a) , and in terms of p(x) and substituting into (4.81a) 
we obtain 


1 c “ - 

- lim zr— L p(t)dtSL (r) [sina (t-x)+ sina (t+x)] 

r^a+0 ^^2 1 

, c (ar) 

- lim /, p(t)dtj Y~? — T [sinct(t-x)- sina(t+x-2b, ) ] da 

r^a-0 “'l 0 i 

= f(x) , bj^<x<c , 


(4.85) 


where the upper and lower signs in the series refer to the symmetric and 
the antisymmetric cases, respectively, and 

, , , K2(ra^)I^(Rg ) - I^(ra^)K^(Rg^) / 

K,(aa )I (Rg ) + I,(ag )K, (Rg ) ° 

Ziiiii Znin 


K, (ra^)I,(Ra ) + I. (ra )K„(Ra^) 
in/n inZn 


K„ (aa„) (Ra^) - I„ (aa^)K (Ra^ 
z n Z n / n z n 


Y for cr^(x) = 0 . (4.86a,b) 


In deriving (4.85) the conditions (4-81b-d) has been used. 


We now observe that for t = x the series and the integral giving the 
kernels in (4.85) are divergent. These divergent parts may be studied and 
separated by considering the asymptotic behavior of the terms in the series 
and of the integrand. In (4.86) note that since R>a, letting r = a+E, for 


large values of a^ we obtain 

K, (aa + ea ) 

1 n n „ -ea. 


= K, (ag ) 
2 ' n 


n 
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(4.87) 


where e Is a small positive number and L^(r)->-l as e->0 and Sim- 

ilarly, note that for large values of a and for a small e- a-r we have 


I^(ar) Ij(aa-ae) 
___ 


IgCoia) 


= e 


-EU 


Thus making use of the following results 


(4.88) 


11m p(t)dtSe ^“slnAn = lim p(t)dt 


3->o 


6-^0 ^1 


sinA 

2(cosh3“COsA) 


= Y p(t)cot dt , 


lim p(t)dt / e ^%inAadot = lim p(t)dt ■ " o' 

£->■0 ° e -^0 ^1 A ^+ e ^ 




(4.89a,b) 


and adding and subtracting the asymptotic values under the summation and 
integral signs in (4.85) we obtain 

r(t+x) 1 { ) 

2b ' 

2b V J 

(4.90) 




1 


t~x t+x-2bi 


1. 


, r ^ Tr(t-x) ' 7T( 

Wi J Cc°t 21, -cot 2b 

2*”^2 I [cosac 2^ + cosac I 


1 rC 


+ — k(x,t)p(t)dt = - y^^fCx) , bj^<x<c 


I, (aa) 


,00 

k(x,t) ~ f (i^~ 77 r's — 1) [sina(t-x) - sina(t+x-2b) ] da 
o \.aaj 


U' lr CO 

+ T— rr- 2[L (a) - 1] [sina (t-x)+sina (t+x) ] 
b^y2 n n n 


(4.91) 


where upper and lower signs and kernels again correspond to symmetric and 
antisymmetric problems, respectively, the kernel k(x,t) given by (4.91) 
is bounded and continuous for all x and t in the closed interval [b^,c], 
and in (4.91), because of uniform convergence, the limit has been put 
under the integral and the summation signs. The integral equation (4.90) 
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must be solved under (4.761), the only condition in (4.76) which has 
not been satisfied. 

At a first glance it may appear that (4.90) is a simple singular 
integral equation of the type discussed in the previous sections. 
However, a closer examination would indicate that the term l/(t+x-2bj^) 
in the kernel of (4.90) is not bounded in the closed interval [bj^,c] 
and becomes unbounded (as 1/x) as x and t go to bj^ together. Hence, 
this term would be expected to influence the singular behavior of the 
solution at the end point b^. Dominant kernels containing, in addition 
to the Cauchy kernel l/(t~x), terms such as l/(t+x~2bj) will be called 
generalized Cauchy kernels . The properties of the solution of singular 
integral equations with generalized Cauchy kernels will be discussed 
in Section 6. 

4.5 Remarks on the Selection of Auxiliary Functions 

In studying multiple series and multiple integral equations if the 
objective is their reduction to singular integral equations, the selec- 
tion of the auxiliary function (i.e., the new unknown function) and the 
procedure followed in the reduction process appear to be quite straight- 
forward. It may be worthwhile to note that in boundary value problems 
in mechanics it is always possible to recognize pairs of "complementary 
functions" bn the boundary having basically the same dimension. Such 
pairs are, for example, the normal and tangential derivatives of the po- 
tential, or the potential and the integral of the normal flux along the 
boundary in potential theory, and the surface tractions and the tangen- 
tial derivatives of the displacements or the integrals of the tractions 
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and the displacements in solid mechanics. If one considers the struc- 
ture of the dominant part of a singular integral equation, namely 

i t^“ ^ (^-92) 

li 

it is clear that (with proper normalizations) the unknown function ^(x) 
and the known function f(x) on the boundary have the same physical dimen- 
sion and in a correct formulation of the problem they invariably are the 
complementary pair on the boundary. Thus, if one does not pay any atten- 
tion to the dimensional consistency in selecting the auxiliary function 
(f>(x) at the beginning, the resulting integral equation may have a sing- 
ular kernel with a singularity either weaker or stronger than the Cauchy 
singularity, l/(t-x). For example, in the former case the kernel is the 
integral of l/(t-x), i.e., log|t-x|, and the integral equation may be 
reduced to the standard form by formally differentiating the both sides 
with respect to x (i.e,, the tangential coordinate), Indicating that the 
particular auxiliary function selected is the complement of f’(x) rather 
than the input function f(x). Analytically, this selection usually does 
not create any difficulty, since one may easily recover and isolate the 
logarithmic kernel by following the procedure described in sections 4.1- 
4.4. However, if the selection is made in such a way that (f)(x) is the 

complement of the integral of f(x), then technically the dominant kernel 

2 

is expected to be the x- derivative of l/(t-x), i.e., l/(t-x) . For a 
Holder-continuous <j)(t) since the integral /4)(t)dt/(t-x) does not exist, 
this formulation becomes meaningless, and besides it is not possible to 
recover the strong singularity l/(t-x) through a normal procedure out- 
lined in the preceding sections. Clearly, the correct thing to do in 
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such a case is to integrate parts of the multiple series (or integral) 
equations so that dimensionally consistent auxiliary function can be 
defined. 

One should again emphasize the importance of writing the boundary 
conditions in limit form on that part of the boundary which will be the 
support of the resulting integral equations. Without this, one may not 
be able to change the order of integrations or integration and summa- 
tion legitimately to evaluate the kernels. Even if this is done, with 
the limit under the integral or summation sign, the resulting infinite 
integrals or series giving the kernels are usually divergent or simply 
meaningless. Consider, for example, the integral equation (4.5) express- 
ing the boundary condition on . If it is not written in limit form 
the kernel becomes 

CO 

K(0,t) = E sin n(t-9) (4.93) 

■ ■ i 

which is not summable. The same thing may be said about the kernels 
arising from the reduction of multiple integral equations or multiple 
series-multiple integral equations (see equations (4.58), (4.85), and 
(4.90)). 

5. NUMERICAL SOLUTION OF SINGULAR INTEGRAL EQUATIONS OF THE FIRST KIND 

In the previous sections it was shown that, unless the problem has 
convective boundary conditions, the mixed boundary value problems in 
mechanics may invariably be reduced to a system of singular integral 

^ ^ ^ Convective boundary conditions generally reduce to integral equations 
of the second kind with a logarithmic dominant kernel which have bounded 
solutions and which may be treated as Fredholm type equations (see Section 
4.2 equations (4,48)-(4.49) and Section 4.3 equations (4.63), and (4.70)). 
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equations of the following general form: 


b. 


Ea..(f).(x) + -^ + k?. (x,t) + (x,t)](|) (t)dt 

I 13 J L ^ 13 J 

= , X = , xeL^ (5.1) 


™i 

where L. = S L... the matrices (a..) and (b..) are nonsingular, the kernels 

1 X3* X3^ X3 

k?.(x,t) consists of terms which become, unbounded as x and t approach 

13 

the end points of and which, with the singular terms b^j/(t-x), con- 
stitute the generalized Cauchy kernels, k^.(x,t) are bounded Fredholm 

X J 

type kernels, and f^(x) are known functions. For the singular equations 
with generalized Cauchy kernels there does not seem to be any general 
method of regularization. The singular behavior of the solution of these 
equations will be studied in Section 6 where a numerical technique for 
solving the integral equations will also be discussed. Also, the treat- 
ment of the singular integral equations of the second kind will be post- 
poned until Section 7. Thus, in this section we vjill consider Only the 
singular integral equations of the first kind with simple Cauchy-type 
singularities which represent by far the largest class of mixed boundary 
value problems in mechanics. The method will be described for a single 
equations defined in the normalized Interval (-1,1), namely 


1 jl Mtl dt + /l'k(x,t)(f)(t)dt - f (x) , -Kx<l . (5.2) 

IT - 1 t-X - 1 

From the development of the method, it will be clear that the extension 
of the method to a system of singular Integral equations (with unknown 
functions (|)^(x) , . ., c|)^(x)) defined in a single interval a<x<b is quite 
straightforward. It is also easy to show that if each equation in the 
system is defined on a (different) union of arcs, the technique developed 
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for (5.2) may easily be used for the numerical solution of this general 
system. For example, consider 


^ k(x,t)(})(t)dt = f(x) , xeL 

li Ij 

n 

L = EL. , L, = (a. ,b.) , b.<a, , 

1 1 ’ 1 ' i’ ’ 1 i+1 


(5.3) 


Defining the following new variables and functions, 


t 2t 

i X 

i b.-a. 

X X 

bi-a^ 

2x 

b.+a. 

i X 

1 D .-a. 

b.-a. 


i i i i 


, a.<t<b. , -Kt.<l , 

X X X 


, a.<x<b, , -Kx.<l , 

’ i X. X ’ 


f(x) = » a^<x<b^ , -Kx^<l , 


k(x,t) = k,.(x,,t.) , a,<x<b. , a.<t<b. , 

xj^ x’ ’ X X ’ j j ’ 


-l<x.<l , -l<t,<l , i = 1, . . ,n 

X 3 ’ ’ 


(5.4) 


and writing (5.3) on each interval xeL. separately we obtain 

b.-a. ^ 

_J 1 


11 

E i /! r 

. . TT -1 b. -a. 
3 = 1 _J 1 




b.+a, b.-a. b.+a. 

t + X - 

2 i 2 2 i 2 


i b.— a. 

+ E /, k. . (x. ,t.)4). (t .) dt. = f.(x.) , 

^ = 13 1 3 3 3 2 3 x^ X 


i = l,..,n , -Kx^<l 


(5.5) 


It may be noted that all the variables x^ and t^ in (5.5) vary between -1 
and 1 hence, the indices i and j 1*® suppressed. Also note 

that in the first terra of (5.5) if i ^ j the kernel is bounded and continuous 
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in the closed interval [-1,1], i.e., for -l<(x^, tj)<l. Thus (5.5) is 
equivalent to the following simple system of singular integral equations 
defined in the normalized interval (-1,1): 

1 ^1 1 

? 4 ~T=r \j(x.t)(}).(t)dt = f.(x) , 

-l<x<l , i = 1, . .n , (5.6) 

where h^^ is the sum of and the corresponding nonsingular kernels in 
the first term of (5.5)^^^. 

Referring to equations (4.10-4.15) and [3] the fundamental function of 
(5.2) may be expressed as 

w(x) = (l+x)“^^^'''°^l(l-x)^^^'’^2 

where and a 2 are (positive, zero, 
is known as the index of the integral equation. The first step in the num- 
erical procedure which will be described in this section is the determin- 
ation of the integers and the index of the problem. As pointed 

out earlier, generally this is not possible without referring to the 
physics of the problem For example, consider the plane contact problem 
shown in Figure 4 for a rigid stamp with a given profile acting on the 
elastic for half plane X2<0, -°°<Xj<=». Using the Fourier transforms, or 
the Green's functions, or the complex potentials [1-4] it can easily be 
shown that in the absence of friction the mixed boundary value problem is 
reduced to the following simple singular integral equation: 

^^^From the analysis and particularly numerical view point another advantage 
of this procedure is that one now is dealing Hith a simple fundamental func- 
tion , for example, of the form w(t) = (l-t2)+l/2 rather than a complicated 
function defined by (4.14) or (4.15). 


, (-l<x<l) (5.7) 

or negative) integers, -(a^ + ~ ^ 
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, -Kx<l 


(5.8) 


i A 4>(t) 

TT -1 t-X 


dt 


f(x) 


where 


2x 


1 b+a 


X = -r -r — , (j)(x) = - a (X|,0) 

b-a b-a yy 1 ’ 


g^vCxj^jO) , a<Xj<b (5.9) 

with Ji and being the elastic constants (K^ = 3-4V for plane strain, 
IC^ = (3-v)/(l+V) for plane stress, V: Poisson’s ratio). Note that the 
fundamental function of (5.8) is given by (5.7) and the solution is of 
the following form [3] 


f(x) = - 


4u 

1+fC 


(|)(x) = F(x)w(x) , -l<x<l 


(5.10) 


where F(x) is bounded in -l<x<l. Thus, the singular behavior of the solu- 
tion is completely determined by that of w(x) . In Figure 4a the contact 
at both ends x^ = a and x^ = b is "smooth" and a and b are unknown. Con- 
sequently at the end points the contact stress <f)(x) must be bounded (and 
necessarily zero). Therefore, in (5.7) we have = 1, = 0, and the 

fundamental function and index become 

w(x) - (1-x ) , K - -1 . (5.11) 


On the other hand, in Figure 4b the leading edges of the stamp are 
"sharp" and a and b are known. Thus, since the contact stress at these 
end points is known to be unbounded, from (5.7) and (5.10) it follows that 

. - 1/2 

= 0 , - -1 , K = 1 , w(x) = (1-x ) . (5.12) 

Similarly for the stamp given in Figure 4c a is known, b us unknown and 
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“l = 0 = a2 


- 1/2 1/2 

K = 0 , w(x) = (1+x) (1-x) 


(5.13) 


Also, for the stamp given in Figure 4d b is knoxm, a is unknown and 

1/2 - 1/2 

= 1 , = -1 , K = 0 , w(x) = (1+x) (1-x) . (5.14) 

Referring now to (4.8) and (4.21) the general solution of (5.8) may 
be expressed as 


(bM = _ w ( x )_ J.1 __f(t)_dt _ 

TT (t-x)w(t) 


+ Cw(x) 


-1<X<1 


(5.15) 


where C is an unknown constant. Although there are very general rules for 
determining the unknown constants such as C, a, and b [3], as pointed out 
earlier in this section, since it is always possible to reduce the general 
singular integral equations to a system defined only in the normalized 
interval (-1,1), for the numerical methods which will be developed in this 
section it is sufficient to state the rules for a simple singular integral 
equation such as (5.8) (or (5.2) for which one simply replaces f(x) by 
[f(x) - /j^k(x,t)(|)(t)dt] ) . One may also note that the statements made here 
for the singular integral equations of the first kind are also valid for 
the equations of the second kind v/ithout any modification. 

(a) K = -1 (Figure 4a) : 

In this case the conditions at infinity require that the constant C 
must be zero and the following consistency condition must be satisfied [3]: 


[f(x) - k(x,t)(^>(t)dt] = 0 . 


(5.16) 


Noting that f(t) contains the constants a and b (see (5,9)), (5.16) provide 
equation for the determination of remaining unknowns a and b. The 


one 
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second equation is obtained by considering the following equilibrium 
condition; 


"•^3 <7yy(xj^,0)dXj^ = c|)(x)dx = P . (5.17) 

(b) K = 1 (Figure 4b) : 

In this case C is the only unknown constant and is determined by sub- 
stituting from (5.15) into the equilibrium condition (5.17) (which must 
be satisfied in all cases) . 

(c) K == 0 (Figure Ac or 4d) : 

In this case again the conditions at infinity require that the constant 
G be zero and (5.15) gives the unique solution without any reference to any 
additional conditions. The problem is solved by assuming that both a and 
b are known and if, instead of the contact area, the resultant load P is 
specified, (5.17) is used to relate the two. 

The numerical methods used for the solution of singular integral equa- 
tions may be considered in two separate categories. The first is a rather 
direct approach which is based on the development Gauss-Jacobi type inte- 
gration formulas for singular integrals. The second is basically a series 
solution with Chebyshev or Jacobi polynomials being/ the related othogonal 
polynomials used in the series expansion. 

5.1 Solution by Gaussian Integration Formulas 

The most common numerical technique to solve a Fredholm-type integral 
equation of the form 

^^^For K = 1 there is always a physical condition such as (5.17) to be sat- 
isfied. For example, in crack problems instead of equilibrium one has the 
singlevaluedness condition (in which P is replaced by zero). 
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k(x, t)(j)(t)dt = f(x) , -l<x<l 


(5.18) 


is the use of some kind of integration formula to evaluate the integral in 
terms of some discrete set of values <j)(tj^) , . , . ,(j)(t^) thereby reducing the 
integral equation to a system of algebraic equations in (f)(t^). In partic- 
ular, if it is appropriate [20,21] one may prefer an integration formula of 
Gaussian type for which (5.18) becomes 

Sk(x^,tj)(|)(tj)W'^ + R^(x^) = f(x^) , i=l,.,,n , (5.19) 

where , (j = l,..,n) is the weighting constant of the related integration 
formula and is the remainder. By selecting n sufficiently large, R^ can 
be made as small as necessary for the desired accuracy and hence, may be 
neglected. In (5. 19) t^, . . , t^ are the roots of the related orthogonal 
polynomial. This highly appealing simple technique could be used for the 
solution of singular integral equations if the Gaussian integration formulas 
for singular integrals were to be available. Some of these formulas will be 
developed in this and the following sections. 

5.1.1. Gaussian Integration Formula for K = 1 

For K = 1 the fundamental function of (5.2) is given by (5.12) which 
is the weight of Chebyshev polynomials (of the first kind) T^(x) . Thus, 
before deriving the integration formula the following property of the Chebyshev 
polynomials will be proved; Let 


T (t. ) 
n k 





r=l,..,n-l . (5.20) 


Then 


n ^^(tJ Jo > j = 0 , 

I ^ = S (5.21) 
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where the Chebyshev polynomials are defined by 


T_^(x) = cos n9 , U^(x) = ■ , cos 0 


n 


X 


(5.22) 


To prove (5.21) note that l'jj(x) and U^(x) are polynomials of degree n and 
consider the following simple fraction expansion: 




n a, 


= S 


1 








Using [22] 


(5.23a,b) 


U 


- T.(x)U„.j(x) - T„(x)Uj.i(x) 


(5.24) 


from (5.23) it follows that 
n(tj^-x) ~ T„(x) 


(5.25) 


lc=l k. ' n 

First part j = 0 of (5.21) follows immediately from (5.20) and (5.25). Sub- 
stituting now from (5.24) into (5.25) We find 


k-i 


U^ (x) - — ^ 


T^(x) 


(5.26) 


which, for x = x^ and 0<j<n is reduced to (5.21) by (5.20), thereby completing 
the proof. 

Let the solution of the singular integral equation nox<7 be of the form 
(5.10) with w(x) as given by (5.12). Let us assume that the unknown bounded 
function F(x) can be approximated to a sufficient degree of accuracy by 


F(x) E i: A,T.(x) 
o J 


(5.27) 
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By using the relation [22] 


1 1 T..(t)dt 


U._j(x) 


j = 0 
. j>0 


(-l<x<l) 


(5.28) 


the singular integral in (5.2) or (5.8) may then be expressed as 


7T -1 t-X 


P 
2 

j=0 


A. i f] 
J 7T -1 


T. (t)dt 




(t-x)(l-t^)^^^ 


S A.U, , (x) 
1 J 3-1 


-1<X<1 


(5.29) 


At X = X , substituting from (5.21) and (5.27), (5.29) becomes 


n E 5 


— J — dt = L L — = ^ ^ 

7T -1 t-X .11 n(t, -X ) . 1,1 n(t, -X ) 

r 3 = 1 k=i k r j=l k=l ^ k r' 




«‘k> 


, n(t, -X ) , , n(t, -X ) 

<^=1 k r k=l k r 


where 


(5.30) 


(2k-l)7T , , Trr , i /r- .n \ 

t, = cos — 7, , k = 1, . . ,n ; X = cos , r = 1, . . ,n-l. (5.31) 

K /n r n 

Note that if the expansion (5.27) is exact, then there is no approxi- 
mation in (5.30) for any n>p. Also note that (5*30) is identical to the 
following standard Gauss-Chebyshev integration formula for bounded func- 
tions [20,21]: 


i — Sd j*. ~ y o(v t ) T (t = 0 

TT ^j_^2yl/2 ^ n ’ n'-V ’ 


(5.32) 


with the important difference that (5.32) is valid for any x whereas (5.30) 
holds only for certain discrete set of x^ given by (5.31) . Using now (5. ID) 
(5,30), and (5.32) the integral equation (5.2) may be reduced to the follow- 
ing system of linear algebraic equations in the unknowns F(tj^) , . . ,I'(t^) : 
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r = 1 , . . ,n-l 


I 


n 

k=i i '■'V'rir + ■■ 

k r 


E ^ F(V = A 
k=l " ^ 


(5.33a,b) 


where A is a known constant and (5.33b) comes from the additional condition 
of the problem (such as (5.17)). 


5.1.2 Gaussian Integration Formula for K - - 1 

For K = - 1 the fundamental function of the integral equation is given 
by (5.11) which is the weight of Chebyshev polynomials of the second kind 
U^(x) . Thus, we will first prove the following property: Let 




(5.34) 


then 




- Tj+iC^r> • 


(5.35) 


Proof: Using [22] 


(l-t^)U^(t) = (n+l)U^_j^(t) - ntU^(t) , 


(5.36) 


the following expansion may be obtained 




(5,37) 


Considering the recursion formulas (5.24) and [22] 


C^(t) = U^(t) - tU^_^(t) , U^(t) = T^(t) + tUj_^j^(t) , 


(5.38) 


from (5. 37) and (5.34) it follows that 


2 

k=l ("+o(tk-^) 


n.““i ““l 


U (x) 
n 


T.^^(xH 


°j(">\+lW 

Un(x) 


(5.39) 


Noting that =0, for x = (5.39) would reduce to (5.35) 

completing the proof. 

Consider now the singular integral in (5.2) with the solution ex- 
pressed as in (5.10) and (5.11). Let the following truncated series 
represent the unknown bounded function F(t) with sufficient accuracy: 


F(t) 2 2 B.U. (t) 
o J J 


(5.40) 


Using (5.35), (5.40), and the relation [22] 


7 4 

7T -1 


u.(t)d-t^)^^^ 


t-x 


dt = - T 


j+1 


-Kx<l 


(5.41) 


for X = X the singular integral may then be expressed as 




i MO . 1 

7T -1 t-X - o 7T -1 t-X 


dt = - S B.T, , , (x ) 
o J J+1 t : 


J-O k=l <"k-\> k=l 


(5.42) 


where 


t,_ = cos ^ = l,..,n;x^ = cos , r == l,..,n+l.(5;43) 


'k 


ni-1 


One may again note that (5.42) is identical to the following Gaussian 
integration formula for a bounded function g(x,t) 

2. 

1 1 o 1 /o ^ 

i g(x,t)(l-t^)'/^dt s 2-^^-^ g(x,t^^ , U^(t^^) = 0 (5.44) 

with the important difference that (5.44) is valid for any x whereas (5.42) 
is valid only for a certain set of x ^ given by (5.43). Also note that 
if (5.40) is exact, again there is no approximation in (5.42) for any n>p. 
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Using now (5.10), (5.11), (5.42) and (5.44) the integral equation 


(5.2) may be reduced to 
2 


n 1-t 
E 

k=l 


n+f ^ ^ ^t" -X ^ = f(x ) , r = i,..,n+l . (5.45) 

k r 


Referring to (5.10), (5.40), and (5.41) and using the orthogonality condition 

fo , k ^ r , 

1 , k = r = 0 , (5.46) 

1/2 , k = r>0 , 




dx 


. ( 


the consistency condition for (5.2) may be expressed as 


[f(x) - k(x,t)(j)(t)dt] 


dx 


1 J.1 dx 


fl dt 


2.1/2 7T -1 .. 2.1/2 -1 t-x 

(1-x ) (1-x ) 


= - /} E B T (x) ^-TJ? = 0 (5.47) 

“1 o 3 3+1 (i-x2)l/2 

Thus, (5.45) implies that the consistency condition of the integral equa- 
tion has been satisfied. Note that in (5.45) there are n unknowns B(t|^) 
and n+1 equations. If K = -1, usually there are two more Unknown con- 
stants such as a and b in Figure 4a and one more condition such as (5.17) 
giving altogether n+2 unknowns and n+2 equations. 

5.1.3 Gaussian Integration Formula for K = 0 

Let the fundamental function and the index of the integral equation 
(5.2) be 

w(x) = (l-t)"^^^(l+t)^^^ , K = 0 . (5.48) 

Expressing again the solution by (5.10) and observing that w(x) is the 

(- 1/2 1 / 2 ) 

weight of Jacobi polynomials ’ (x) , it will be assumed that the 

unknown bounded function F(t) may be approximated to a sufficient degree 
of accuracy by 
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(t) , -Kt<l . 


(5.49) 


F(t) E 

“ o 3 J 

From (5.10), (5.49), (5.48), and using the following general property of 
Jacobi polynomials [23] 

i ;l p<“.6)(t)(l-t)“(l+t)® ^ 2-V:“'-6>W , 


•Kx<l , K = -(a+8) = (-1,0, or 1) , -l<(a,3)<l , 


(5.50) 


the singular integral in (5.2) may be expressed as 


i /J dt 

TT -1 t-X 


o J 3 


-l<x<l 


(5.51) 


Let 


F^( = 0 , k = l,..,n , 


(5.52) 


and consider the expansion 

n 3 n 3 


n c. 


= E 


, 3<n . 


1 ti.-x 

Using (5.24) and the relations [21] 


(5.53) 


(-1/2, 1/2) . ^^2_ . ^ r(n-fl/2) 1 

n ^1/2 z ^2n+l^^^ V 


P ^^^^’“^^^^(23^-1) = u (z) , 

n 2n 


(5.54) 


equation (5.53) may be modified as 
Uo„_o._i(y) " c. 


r(.i-H/2) ^2n-23-l^"^ _ “ 

•" * t ^ ZI / "\ ' ' " *-> 


To .,(y) 1 t,-x 


, X - 2yM ■ 


j: TT -■ -2n+r^/ ^ nc 

Indicating that the expansion in (5.53) is indeed possible. 
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(5.55) 



Referring now to the results in Section 5.1.1 and to (5.54) and using [23] 


d -p (a, 3). V n-t-g (a+1,3-1). v _ 

dt n 1+t, . n r 

k 


k 


= 0 » k - l,..,n, (5.56) 


the coefficients c, are found to be 
k 

, 2(14-t^) (.1/2, 1/2) 

^k 2n+l 3 ^ k'' • 


(5. 57) 


If we let 


P (l/2»-l/2) y ^ 0 ^ 

n ■ r 


r = 1, . . ,n 


(5.58) 


from (5.53) and (5.57) it follows that 


n 2(l+t, ) 
S 

k=l 


(-1/2.1/2), , . (5.59) 


2n+l 3 


Finally, from (5.49), (5.51) and (5.59) we find 




7T -1 t-X TT -1 t-x 1-t . 2n+l 

r r K-i 


t, -X 

k r 


(5.60) 


V7here 


(- 1 / 2 , 1 / 2 ) 

n ''^k^ 


2k- 1 

0 , t.j^ — cos ( 2 .J 2 -}-'! ’ k — 1, . . ,n 


n r 


, 2r7T . , 


(5.61) 


Equation (5.60) too is identical to the corresponding GausS-Jacobi integra- 
tion formula fot a bounded function g(x,t) given by [21] 


1 1 14 .^ 1/2 n 2(t, +1) 

^7 8(x,t)(i±|) dtsE 

k=l 


(5.62) 
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with again the important difference that (5.62) is valid for any x 
whereas (5.60) is valid only for x = x^, r « l,..,n. 

Using the integration formulas (5.60) and (5.62) » the singular in- 
tegral equation (5.2) may easily be reduced to the following system of 
linear algebraic equations in F(tj^) , . . ,F(t^) ; 
n 2 ( 1+t, ) . 

f , 2n+l ° (5.63) 

k=l k r 

where t, and x are given by (5.61). 

xC IT 

1/2 1/2 

If the weight function is w(x) = (1-t) /(1+t) following a simi- 

lar procedure, it can be shown that 


1 .1 FCO A-U 

TT -1 t-X ''1+t^ 


1/2 


n 

dt = S 


k=l 


2n+l t, -X 
k r 


(5.64) 


and (5.2) is reduced to 


n 2 (1-t.) , 

k=l k r 


r = 1, . . ,n (5.65) 


where 


p (1/2, -1/2). . 

n '‘^k- 


0 , tj^ = cos(-~y) , k = 1,. . ,n 


P ^ 1/2, 1/2)^^ ^ ^ ^ 

n r^ ’ 


/2r-l „v .1 


(5.66) 


5.2 Solution by Orthogonal Polynomials 

In the previous sections it was indicated that the fundamental function 
of the singular integral equation (5.2) 1$ of the general form 

w(x) = (l-x)“(l+x)^ , K - (o+B) = (0,+l) . (5.67) 

Noting that (5.67) is the weight of Jacobi polynomials P^ ’ (x) , it is 
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natural to look for a solution of (5.2) in the following series form: 


{j)(t) = w(t)S C P , -Kt<l . 

o ^ 


(5.68) 


Thus, substituting from (5.68) and (5.50) into (5.2) we obtain 


00 .-K (-a, -3) 

2 C [- P 


n=0 


n sxrnra n-K 


(x) + k^(x)l = f(x) , -Kx<l , 


k^(x) = k(x,t)P^^°^’^\t)w(t)dt . 


(5.69) 


Equation (5.69) can further be reduced to an infinite algebraic sys- 


tem in C as follows : 
n 

o-K 


" sin-rra ^ . ^kj^j “ % ’ O.l,-- 

3=1 

where for K = -l, C_j^ = 0 and 

^kj -1 ^^(x)k^(x) (l-x) “‘(l+x) '^dx , 

c^ " f(x)P^^ ^\x)(l-x) '^(1+x) ^dx . 


(5.70) 


(5.71) 


the constants 0^^ come from the orthogonality condition 


p <^“»^^(t)P, ^°^’^\t)w(t)dt 

JL 1TL IC • 


= /“ ’ " 


^ k 


e^^(a,3) , n = k , 


k = 0,1,2 




0„(a,3) 


e^(a,3) 


[\ w(t)dt .■2"®^’rca«)r(B±iI 


r (ot+3+2) 


2 °’~^^'i'ir (k-fct+i) r (k+g+i) 
( 2 k+a+ 3 +i ) kl r (k+a-l- 0 + 1 ) 


k = 1,2,. 


(5.72) 


The infinite system (5.70) may be solved by the method of reduction [24]. 
Note that if K = 1 in the reduction letting k = 0,.. ,n in (5.70) involves 
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the unknovms , . . , , that is , there is one more unknown than the 


number of equations. The additional equation is provided by the follow- 
ing physical condition (e.g., (5.17)): 


.1 


/j c}>(x)dx = A = C^9^(a,3) 


(5.73) 


where A is a known constant. On the other hand if K = - 1, from (5.69) 
and (5.71) it is clear that the first equation in (5.70), (i.e., k = 0) is 
equivalent to the consistency condition (5.16). It should further be noted 
that the integrations necessary for the evaluation of the constants Cj^^ and 
Cj^ in (5.71) are of Gauss-Jacobi type may easily be evaluated by using [20] 

g(x)(i-x) °^(l+x) ^dx = S \s(x^^) . 

k=l 

P (““’"^)(x, ) = 0 , k=l,..,n , 

n k 


w _ (2n-a-3+2)r(n-a+l)r(n-3-H) * 

\ “ (n+l)!(n-^-3+l)r(n~a-3+l) 


,-a-3 


4 - P ^"“’"^^(x^) 


dx n 


V n+1 


(5.74) 


2 “ 1/2 

In the special case of K = 1, w(x) = (1-x ) the related ortho- 
gonal polynomials are T^(x), (n = 0,1,2,..), and we have 

CO . ' ■ . ■ ■ 


(f)(x) = (1-x^) S A T (x) 
o 


^ ll U^(t)U^^(t) (l~t^)^-^^dt =' 


0 , n k 

tt/P , n = k 


(5.75) 

(5.76) 


TT 

2 


n=0 


(5.77) 
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(5.78) 


f (x)Uj^(x) (l-x^)^^^dx , 

\n " -1 \(x)(l-x^)^^^dx /J k(x,t)T^(t)(l-t^)"^^^dt 
(j)(x)dx = A = ttA^ 

where (5.28) and (5.46) have been used. Truncating (5.77) at the Nth 

term and retaining first N equations, (5.77) and (5.78) give the N+l 

unknown constants A ,A, , . . ,A . 

o 1 ’ n 


In the other special case K= - 1, w(x) 
be shown that 




it may easily 


^(x) = (1-x^)^^^ S B U (x) , 

o n n 


- T \-l + S-O \n\ = \ • 

\ = l\ f (x)Tj^(x) (1-x^) , 

b^^^ = T^(x)(l-x^) ^^^dx k(x,t)U^(t) (l-t^)^^^dt 


(5.79) 

(5.80) 


n 


where (5.41) and (5.46) have been used. Again in (5.80) the first equa- 
tion (k = 0) corresponds to the consistency condition (5.16). 


Some typical applications of the numerical methods described in this 
section may be found in [25-27] which also Include extensive references to 
the solution of mixed boundary value problems in mechanics obtained by 
applying these methods. 
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6. INTEGRAL EQUATIONS WITH GENERALIZED CAUCHY KERNELS 


As pointed out earlier in this article, in some mixed boundary value 
problems formulated in terms of singular integral equations, in addition 
to Cauchy type singularities the kernel may contain terms which become un- 
bounded as both variables in the kernel approach an end point. For example, 
in the torsion problem shown in Figure 3 and formulated by the integral 
equation (4,90) if b>c, separating all the singular terms, (4.90) may be 
expressed as 

? ~ t+x-2bj^^ 


+ kj^^(x,t)p(t)dt = - ]ij^f(x) , bj^<x<c , ’C6I1) 

where kj^^(x,t) is a Fredholm type kernel which is bounded in the closed 
interval bj^<(x,t)<c and is the sum of k(x,t) (see, 4,91) and the nonsingular 
terms in cotangent or cosecant kerviels appearing in (4.90). Note that in 
(6.1) the term l/(t+x-2bj^) becomes unbounded at the end point bj^. 

Similarly, if c = b the kernel in (4.90) contains additional "singular” 
terms which may be separated by observing that 


JL 

2b 


cot 


'TT(t-x) 

2b 




JL cot + h (x t) 

2b 2b t+x-2b • ’ 


2b 2b 


Tr(t-x) 1 




IT F(t4-x) 

2b 2b 




(6.2a-d) 


where are bounded in the closed interval bj^<(x,t)<b. Thus, sub- 

stituting from (6.2) into (4.90) we find 


lb ^*11 

-/ p(t)dt[(,i +-) - 


-X t+x-2bj^ U 2 t+x-2b 


+ k2j(x,t)p(t)dt = - 


(6.3) 


where, again is a bounded kernel. As seen from (6.1) and (6.3), the 
integral equations for the symmetric and the antisymmetric problems have 
the same dominant parts but different Fredholm kernels. In (6.3) in addi- 
tion to the Cauchy singularity, the dominant part of the kernel contains 
terms which become unbounded at both ends, b^ and b. The dominant kernels 
of this type are defined as generalized Cauchy kernels . 

In the torsion problem under consideration there is one more limiting 
case which is worth considering. This is the case of bj^ - 0, (e.g., a 
broken shaft). Thus, for bj^ = 0, c = b, observing that in addition to 
(6.2a and c) , around the end points one may write 


F _ TT(t+x) 1 , 1 . N 

2b 2b t+x t+x-2b ’ 


TT 7r(t+x) 1 1 

2b 2b t+x t+x-2b 


+ hg(Xjt) , 


(6.4) 


the integral equation (4.90) may be expressed as 


lb ^1 1 ^11 


+/ k (x,t)p(t)dt = - y,f(x) , 0<x<b 

0 Jt i 


(6.5) 


where h^, h^ , and k_^ are bounded in 0<(x,t)<b, and the upper and lower 

D O JX ~ 

signs in (6.5) correspond to the symmetric and the antisymmetric problems , 
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respectively. Structurally, (6.5) Is identical to (6.3), however, at 
X = 0 their solutions may have quite different behavior. 

6.1 A Plane Elasticity Problem for Nonhomogeneous Media 

Perhaps the most typical application of the singular integral equa- 
tions with a generalized Cauchy kernel arises in the study of crack prob- 
elms in nonhomogeneous elastic solids. Consider, for example, the three- 
dimensional problem in which part of periphery of a plane crack extends to 
the bimaterial interftice in a nonhomogeneous medium. Consider the cross- 
section shown in Figure 5a. In the neighborhood of the point of interest 
0 let the interface be a smooth surface and, for simplicity, let the crack 
plane be perpendicular to the interface. Through a proper superposition 
one can separate the singular or the perturbation part of the problem in 
which statically self-equilibrating crack surface tractions are the only 
external loads. From the viewpoint of the stress state around the point 0, 
the perturbation problem in turn may be assumed as having three components, 
namely, the in-plane extension (mode I), the in-plane shear (mode II), and 
the anti-plane shear (mode III) with the corresponding surface tractions 
shown in Figure 5a. Clearly, the behavior of the solution of the anti-plane 
component of this problem around 0 will be identical to that of the torsion 
problem shoTTO in Figure 3 at x = 0 = b^^, r = a. This behavior should be 
completely characterized by the first two terms of the genralized Cauchy 
kernel appearing in (6.5). 

Similarly, the singular behavior of the solution under ip-plane loading 
conditions will be equivalent to that of the plane strain problem for two 
semi- inf inite media having a finite crack perpendicular to and ending at the 
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interface shown in Figure 5b. Formulation of this problem is straight- 
forward. For example, referring to [28] for details, in the case of in- 
plane tension problem, i.e., for the following symmetric loading conditions: 


a^ggCr,?!) = f(r) , ^ ’ 0<r<b , 


the integral equation of the problem may be expressed as 

1 r 1 . ""l . V . ""3^ V 


7T o '■t-r t+r . , .2 .3 

(t+r) (t+r) 


-](})(t)dt = f(r) 


2y 


1 


0<r<b , 


where 


4>(r) = - [ujQ(r,TT+0) - u^ (r,TT-0)] , 

Cf =Y~ mCl+Kj)/ [2(m+l<2)] - 3(l-m)/ [2(l+tnKj^)] , 

=-- 6(l-m)/(l+tniCj^) , == 4 (m-1) / (1+mK^) , m = 


( 6 . 6 ) 


(6.7) 


( 6 . 8 ) 


and and are the elastic constants (i = J,2, = ^~^i plaae strain, 

= (3-V^)/(l+v^) for plane stress, Poisson’s ratio). For the in-plane 
shear problem one obtains the same (generalized Cauchy) kernel as that shown 
in (6.7). Thus in the general problem of more complex geometry and loading 
conditions the system of three integral equations representing the problem 
will have the dominant kernels given in (6.5) and (6.7) (in uncoupled form) 
and coupling will be through Fredholm type kernels only. 


Referring to Figure 5 if we now assume that the crack crosses the 
boundary and extends into the medium 2 (dashed lines) , it is again not diffi- 
cult to argue that the dominant parts of the related system of integral equa- 
tions will be the same as that of the idealized through crack problem for 


EBPRODUCIBILffY OF THV 

OEffilNAL IS PO'M 


two elastic half planes shown In Figure 5b. The derivation of integral 
equations is again straightforward. For example, in the symmetric problem 
selecting the unknotm functions 

g 

^l(r) ==-37 [uj^g(r,7T+0) - u^g(r,TT-0)] , 0<r<b^ , 



+0) - U2g(r,-0)] , 

0<r<b2 , 

(6.9) 

the following crack 

. surface tractions 



CieeCrj-fr) = fj^(r) 

’ CTirQ(r,Tr) = 0 , 

0<r<b j , 


<7200(^>O) = £3(1^) 


0<r<b2 , 

(6.10) 


the integral equations for the perturbation problem may be expressed as 
(see [29] for details) 


1 b. <{>-,• (t) 2 

i. f X __i 


. bj 1+K. 

/ 1 dt + E k (r,t)cj) (t)dt = -x— f . (r) , 

TT o t-r j=l ^ o 13 X 


i = 1,2; 0<r<b , , 

X 


( 6 . 11 ) 


where 


k-1 

3 c.,r 

k. (r,t) = E — r , i = 1,2 , 

k=l (t+r) 

j k-1 

2 d., r 

k..(r,t) = E ^ , i,j = 1 , 2 ; i 7^- 3 

k=l (t+r)*^ 


(6.12a,b) 


c^j^ = 1/2 - m^(l+ic^)/E2(m^+lC2)] - 3(l-m^)/[2(l+mj^ic^)] , 

C 12 = 6 (l-mj)/(l-hn^K^) , Cj^^ = 4(nij^-l) / (l+m^l^k^) , 


C01 = 1/2 “ m^(l+K^)7[2(m2+K^)] - 3(l-m,) / [2 (l+m,l<c^) ] , 


*21 


2 2' 
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^11 6(l-m2)/(l+ni2K2) » <^23 ~ ’ 

^11 ~ 3(l+'^j)/[2(m2+Kj)] - (1+Kj^) / [2(l-hn2K2) 1 » 

^12 ~ (^-+Kj^)/( 1+1112 1 ^ 2 ) ~ ( 1+1^2^ )/(ni2'*^'^l^ ’ 

^21 " 3(l+K2)/[2(mj^+K:2)i ~ (I+K 2 )/ I^^Cl+irij^Kj^) 1 , 

^^22 ~ (I+IC 2 ) - (l+lC2)/(mj^+K2) , 

= ^2''^! ’ ™2 "" ^1^^2 * " (6,13) 

From (6.11) and (6.12) it is seen that the integral equations contain dom- 
inant kernels only which are of generalized Cauchy type. 


6.2 The Fundamental Functions 

The integral equations (6.5), (6,7), and (6.11) are some special cases 
of the following system of singular integral equations with generalized 
Cauchy kernels: 


1 N b A k jk -1 

Ly f ^ y B (x_a ) (t-z. ) ^ 

TT n=l a„ t-x ic=0 mnk n , k In 
n dx 


*^n 3 H ^ ~ ^ 

+ 2 C . (b -x) r (t-z„ ) ](|) (t)dt 

mnj' n ^^3 2 n" ^n" 


N b 


n 


n=l 


where 


z , = a 

In n 


^^(t)dt 

= 

m 

, m = 1 , . . ,11 , a <x<b , 

’ ’ ’ ’ m m ’ 

iSn 

o<e < 2 f 

9 


n 

■ ^ j 

-F<W <1T 

Tl 

) Tl— j 


^2n n ' " n 

A , B ,, C . are known constants, f. ,..,f„ are known input functions, 
mn mnk mnj i N 


(6.15) 


k (x,t) are Fredholm kernels, and 6, are the unknown functions, 

mn 1 

In most practical problems 9^ = 'H' and 0)^ = 0. However, occasionally one 

may encounter cases in which 6 = + Ti/2, w = + tt/ 2 with the related terms 

n n 

in the kernels of the form [28] 

(t-a )/[(t-a + (x-a )^] , ' (b -t)/ [(b -t)^ + (b -x)^l . (6.16) 

n n . n n n n 

Referring to Figure 6a and (6.15) note that as the variable x varies on the 
line of integration L^, the variables and vary on lines and 
respectively. In some problems the kernels may contain terms of the form 
(6.16) as vjell as, for example, (6.12). This means that in Figure 6a at a 
given end point there may be more than one auxiliary line or i2n’ 
meaning that there may be additional terms in (6.14) and (6.15) defined by 
additional angles 0 or w . 

In this type of problems another point which requires special emphasis 

is whether the cuts L = (a ,b ) , (n = 1,..,N) intersect each other or not. 

n n n 

The only type of intersection which is physically meaningful is for two 
adjacent cuts to have a common end point (Figure 6b). A physical example 
for this would be the through crack problem shown in Figure 5b. In such 
problems, aside from the necessary changes in the limits of integrations, 
the integral equations (6.14) remain unchanged. However, as will be pointed 
out in this section, one may have to be somex^hat more careful in determining 
the fundamental functions. In the following, first it will be assumed that 
all end points are distinct, then an example for coinciding ends will be 
considered. 

To find the fundamental functions of (6.14) we let 
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ct 3 

= 8„(t)w (t) = g^(t)(b„ -t) ’^(t-a^) ” 


n 


n n 


"n n 


n 


- l<Re(a^,3^)<0 , n=l,..,N , g^(a^) 0 , ^ ° 

where the function g (t) , (n = 1,..,N) is Holder-continuous in a <t<b , 

n n- - n 

GLy S 

(t-b^)^(t-a^) ri is any definite branch which varies continuously in 

a <t<b , and a and 3 are unknox'jn constants which have to be determined, 
n n n n 

Consider now the following sectionally holomorphic functions: 


g„(c) 


n 


^n 


n 


z (6.18) 


Examining the singular behavior of F^(z) around the end points (see, e.g., 
[3] , Chapter 4) and separating the principal parts we obtain 

T^/\ / \/i \ 0 II / « T1 

F (z) = - g (a )(b -a ) “ — : ■■ ■ w (z-a ) 
n ^ n n • sinTry 

o a 

3 1 a 

+ g (3 ) (b -a. ) ^ — — -Q— (z-b ) ’^ + G (z) , 

n n n n sxmr3 n n 


n=-l,..,N . 


(6.19) 


The functions ^^^(z) is bounded everyv/here except possibly at the ends 
near which it may have the following behavior: 


G^(a) < 


n 


z-a 

' n' 


— , p < - Re (3 ) , G (z) < 

Pr, n n^ ’ ' n ' 


n 


z-b I ^n 


n' 


r <-Re(a ) 
n n 


( 6 . 20 ) 


where A , p , B , and r are real constants, 
n ^n’ n’ n 


Using the Plemelj formulas (3. 61), from (6.19) we find 


1 b 14. 

dt=ftF/« + F„-wi 
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( 6 . 21 ) 


CL 0 

g„(a )(b -a_) ”cot(Tr6 ) (x-a ) ^ 


■"n ' n ' n n 


n n 


3. a 

+ 8„(b )(b -a ) "cot(7ra ) (b^-x) ” + H (x) , 


n n n n 


n' n 


n 


n = 1 , . . ,N , a <x<b , 

where H^(x) is bounded in a^<x<b^, and at the end points its behavior is 

similar to that of which is given by (6.20). Prom the definition of 

(6.18), and from Figure 6a it is clear that for z = z, ^ a and 
n In n 

z = ^ 2 ^ i.e., on and ^ 2 ^, f^(z) is holomorphic. Therefore 

1 b 

7T -^a.^ t^^^z " ^n^^sn^ ’ ^sn^^sn ’ s = 1 ,2, n = 1 , . . ,N , (6.22) 

n sn 

or, substituting from (6.19), using (6.15) and separating the principal parts 


1 b <i> (t) a ^^n^®n 

_ _/• n _B dt = - g (a ) (b -a ) ^ — — ■: — 5 (x-a ) ^ 

ir a_. t-Zi _ °n n n n sinTTp- n 


3. 


A u. ^ 

n In 


n 


+ h (x) , a <x<b 

n n n 


1 V 

dt 

F a t-z„ 
u 2n 


n iW a 

3 n n a 

g (b ) (b -a ) ^ — (b -x) + s (x) , 

n n n n sxnFO n n 

n 


a <x<b 
n n 


(6.23) 


where h (x) is bounded in a <x<b , s (x) is bounded in a <x<b , and the 
n n ~ n n n- n 

behavior of h and s near the ends a and b , respectively, is determined 
n n n n ’ 

by that of G^(z) as given by (6.20). 

Since f’^^Cz) is holomorphic on and l 2 ^ (see (6.18) and Figure 6a) , 
using (6,23) we may also write 


— / (t)(x-a ) — 7- (t-z. ) dt = (x-a ) — r- F (z, ) 

F a ^n^ n , k In n , k n In 

n dx dx 
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- - S„(%)<V%> " ^ sintf 3„(6„-l)---(g„-k+l)(x-a^) " 

n 


(x 1 , (x) , k 1» • • jK , a <x<b_ , 

“ ti n n n 


i / " ♦ (t)(b -x)J -ir Ct-a,„)-ldt - (b -x)j ^ F_(x,_) 


7T a ’^n' n 
n 


j ^n^‘'2n^ 


o iw a 

. 3 n n ( 

= (-l)^g (b )(b -a ) ■ —■• a (ot -1) • • • (a -j+l)(b -x) 

°n n n n sinTTa n n ^ n 

n 


+ (b -x)J A 


+ (b^-x) — j s^(x) , j = , a^<x<b^ , n = l,..,N . 

cL^c 

(6.24) 

Equations (6.21), (6.23) and (6.24) give all the terms in the dominant 

part of the integral equations (6.14). Thus, by substituting from (6.21), 

(6.23) and (6.24) into (6.14), observing that for the rath equation a <x<b , 

m m 

multiplying both sides first by (x-a )~^"i and letting x->-a , then by (b -x)~° 

m ° m m 

and letting x>b^, and also observing that g(a^) 7^ 0, g(hj^) ^ 0, 0<r^<Re(-a^) 

0<p <Re(-3 ), we obtain 
~ n n 

13 (0 "7t) Ktm 

A cos7t3 + e “ [B +Z B ,3 (3 -1) * ’ (3 -k+1)] 

mm m mmo mmk m m m 


= 0 , m = 1,..,N , 

i 

A cosTra + e [C +E C . (-l)'^oi (a -1) • • * (a -k+1)] 

mm m mrao -js! mmj m m m 


=0 , m=l,..,N , (6.25a,b) 


The characteristic equations (6,25) provide two sets of (highly nonlinear) 


algebraic equations to determine the unknoton constants O.^ and 3j^» (m - 1,.., 
With ot and 3 determined, (6.17) gives the fundamental functions w (t) , 

(m = 1, . . ,N) . 
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In the foregoing analysis it was assumed that the 2N end points a^, 
are distinct. In practice this may not always be the case. To give 
an example in which some of the end points may coincide, consider the 
case of, for example, in (6.14) (Figure 6b). As seen from Figure 6b 

(see also Figure 5), physically the problem really has 2N-1 points of sing- 
ularity, and mathematically there are only 2N-1 irregular points. This 
implies that in the expressions of the fundamental functions there will be 
only 2N-1 exponents, i.e., Hence, the fundamental functions of 

(6.14) given by (6.17) will have to be modified as 


W^(t) 

“1 ^1 

= (bj^-t) 

, a^<t<b^ , 


WjCt) 

a 

= (hp-t) (t-bp 

, b^<t<b2 , 


w„(t) 

a 6 

= (b^-t) ""(t-a^) ^ 
n n 

, , n = 3, . ,,N 

(6.26) 


The solutions of (6.14) are still of the form d) (t) = g (t)w (t) , and the 

n n n 

conditions on ot , 3 , g (a ) , and g (b ) , (n = 1,..,N) shown in (6.17) are 
n’ n’ ^n n ’ ‘’n n ’ > > / 

still valid. AlsjO with a„ = b, and So = Ot, the definition of F (z) , (6.18) 

z i Z i n 

and the asymptotic expressions (6.21), (6.23), and (6.24) are still valid. 
Thus, proceeding as before, one obtains (6.25a) for m = 1,354,..,F[ and (6.25b) 
for m = 2,3..,N giving altogether 2N-2 algebraic equations for 2N-2 unknowns 
^1 ’ ^3 ’ * * ’ and 0 I 2 > ‘ * 


Writing (6.14) for m = 1 and letting xr^bj^ and for m = 2 and letting 
x^ap = bj^ one obtains two homogeneous linear algebraic equations for gj^(b^) 


and g 2 (bj^) of the form 


S Ci^(a^)s,(bp - 0 


1.2 


(6,27) 
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Since g^(bj^) ^ 0 (6.27) gives the following characteristic equation to 
determine : 

|c_(apl = 0 . (6.28) 

Note that gj^(bj^) and g 2 (bj^) are not independent and are related through 
(6.27). This provides the additional condition which was eliminated by 

( 1 ) 

letting = bj^ and which is necessary for the unique solution of (6.14) 


As an example for coinciding end points considet the crack problem 
formulated by (6.11) and (6,12) and shown by Figure 5a (with dashed lines 
included). Note that this problem was formulated by using polar coordinates. 
Hence in the terminology of this section, a^^ = 0 = 82 j 3^^ = 32 "^3, and the 
fundamental functions and the solution may be expressed as 

a o -TTia a ^ 

0<t<b , -l<Re(a ,3)<0 , g (b ) ?^ 0 , g (0) ^ 0 , n = 1,2. (6.29) 

n n n n n 


Consider now the following sectionally holomorphic functions 

-TTia 


Fn(.) 


_1 

TT ■ O 


/ ^ -J dt = / ’^(t-b ) t*^ -5 dt 

t-z TT o n t-z 


n =* 1,2 


(6.30) 


where, noting that ( 6 . 11 ) is derived in polar coordinates, the complex varia- 
ble z = r+ip for n — 1 and n = 2 is defined in such a way that in each case 
the cut (0<r<b^) lies along the positive real axis. Repeating the analysis 


^ ^For example, in collinear crack problems letting a 2 =b]^ essentially elimin- 
ates one crack while the number N of integral equations and unknown functions 
remain the same. Physically, since there are N-1 cuts, there will be only 
N-1 single-valuedness conditions. However, the general solution of N integral 
equations will have N arbitrary constants. Thus for a unique solution one 
more condition is needed which is provided by (6.27) (see [29]). 
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given previously, we obtain 


ct -fri 3 p, ft 1 

F (z) = - g (0)b ” z*^ + g (b )b‘^ (z-b ) + G (z) , 

n n sxmr3 '^n n n sinira n n 

n 


a 


a 


, b (j) (t) V. o _ 

— f — ^ dt = - g (0)b "cotTTgr^ + g (b )b^ cotTTa (b -r) “ + H (r) 

ir o t-r ^n n ®n n^ n n' n n 


1 b 

TT 0 t+r 


a 


± f n _ii_ ^ p (_r) 


n 


- g(0)b ^ •■ — . -—- Q + u (r) 

. n sinir 3 on 


^ ° (t+r)^ 

^ f ^ . - o - dt 
"" ° (t+r) 3 


r 4 - F,(-r) = g^( 0 )b^ " r® + u„(r) , 


dr n 


‘In' 


2 ,2 

— ~ — F (-r) = 

2 , 2 n'' '' 

dr 


a 


“2n« • 


0<r<b , n = 1 , 2 , 

n 


(6.31) 


where the behavior of H (r) and u. (r) , (n = 1,2; j = 0,1,2) around r = 0 

n jn 

is similar to that of G (z) around a = 0 as given by (6.20), otherwise they 

n n 

are bounded functions. 

If we substitute from (6.31) into (6.11), multiply both sides first by 

_ct “6 

(b -r) n and let r+b , then by r ^ and let r->-0 we obtain 
n n 


cofrra = 0 , n = 1 , 2 , 

n 


a 

[cosTTg + - 3 c^ 2 + 3(3-1) .;^/2] [g^(0)bj^ /sin'n'31 


+ (d^j - 3dj_2) [g2(0)b2 /sinF3] = 0 , 

a, 

(di2 - 3d22)tgi(0)bj^^ Vsin7r3] 




(6.32) 


+ [cos 7 r 3 + ^2^ - 322 + 3(3-1)022/21 [g2(3)b2 /sin7T3] - 0 . ( 6 . 33 ) 
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Equations (6.32) are the well-known characteristic equations for crack 
tip singularities which are fully imbedded in a homogeneous medium and 
give “ 0^2 “ -1/2. Writing the determinant of the linear system (6,33) 
one obtains the third characteristic equation giving 3 as follows: 

[cosTTp + 3(3-1)0^2/21 [costt 3 + C2J^ - 3C22 

+3(3-1)023/2] - (djj-3d^2)(d2i-3d22) =0 . (6.34) 

It should be noted that as expected (6.34) is identical to the character- 
istic equation giving the stress singularity at the apex of two bonded 
quarter planes by using an entirely different method (see, e.g., [30]). 

Going now back to the other examples discussed in this section, in the 
torsion problem for c<b the dominant part of the integral equation is given 
by (6.1) which is a very special case of (6,14). Thus defining the solution 
and the fundamental function by 

p(x) = g(x)w(x) , w(x) = (cKx)“(x-bj^)^ , bj<x<c , (6.35) 

and applying the procedure described in this section we find 

cos'rra =0 , cos7r3 = • (6,36) 

If c = b, there is a contributiQn to the genera3.ized Cauchy kernel at both 
ends, the integral equation is given by (6.3), and characteristic equations 
are found to be^^^ 

cosTO = yj^/(y^+iJ2) , 003773 = 112/(1^1+112) • ( 6 . 37 ) 

In the case of the "broken shaft" that is, for b^^ = 0, c = b, the integral 

(^) In this case note that in (6. 14) N=lj Kj^=0=Jj^, 0j^=7r, u)j^=0. 
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equation is given by (6.5) from which the characteristic equations are 
obtained as 

cosTra = , cosTTg = (y 2 +yp/(y^+U 2 ) » (6.38a,b) 

where the upper sign is for the symmetric and lower sign is for the anti- 
symmetric loading. Note that for symmetric loading (for which cjQ^(a,x) = 
p(x) = p(-x). Figure 3) 3=0, that is the stress state at (r = a, x = 0) 
is bounded, whereas for antisymmetric loading the characteristic equation 
(6.38b) is identical to that of a semi-infinite crack perpendicular to a 
bimaterial Interface under anti-plane shear loading. 

For the plane problem of a crack perpendicular to and terminating at 
the interface of two bonded half planes (Figure 5) , the integral equation 
is given by (6.7). Thus, by defining, 

<l>(t) = g(t)w(t) , w(t) = (b-t)\^ , (6.39) 

the characteristic equations are found to be 

cosm = 0 , 2 Yj^cos 7T(3+1) + Y2(3+1)^ “ ^3 ^ (6.40a,b) 

Yj = (m+K^) (1+mKj^) , Y2 - m = y2/^l ’ 

Y3 = (l^m) (mfK^) + (1+mK^) (m+Kp) - ra(l+ic^) (I+K^^) « 

Again (6.40a) is the well-kno^«m result giving a = -1/2 and (6.40b) is identical 
to the characteristic equation for a semi-infinite crack perpendicular to and 
terminating at a bimaterial interface under in-plane loading (see, for example, 
[28]). 
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6.3 Numerical Method for Solving the Integral Equations 
with Generalized Cauchy Kernels 

As shown in Section 5 of this article, a general system of singular 
integral equations such as (6.14) can always be expressed by means of a 
simple system in which both variables x and t vary in the normalized in- 
terval (-1,1). Furthermore it was also indicated that, once the method 
of solution is developed for a single eqtiation, it can easily be extended 
to a system consisting of any number of equations. The numerical method 
for solving the singular integral equations with generalized Cauchy ker- 
nels will therefore be described for a single equation of the following 
form only: 

~ + k^(x,t)]c|)(t)dt = f(x) , -Kx<l , (6.41) 

— 1 * 

where f is a known function, k„ is a kno^m Fredholm kernel, k (x,t) becomes 

c s 

unbounded as x and t appraoch the end points + 1, is otherwise bounded, and 
with 1/ (t-x) form the generalized Cauchy kernel. Generally, (6.41) must be 
solved under an additional (physical, such as an equilibrium or a single- 
valuedness) condition of the form 

^(t)dt -A , (6.42) 

where A is a known constant. As indicated previously, the unknown and the 
fundamental functions may be expressed as 

(f)(t) = g(t)w(t) , w(t) - (l-t)°^(l+t)^ (6.43) 

where C'- and ^ are. known constants with -l<Re(a,3)<0 and g(t) is an unknown 
function which is bounded in Observing that w(t) is the weight 

function of Jacobi polynomials P^^^’^^(t), the integral equation may be 
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solved by using a numerical method based on a Gauss-Jacobi integration 
formula which is similar to the methods described in Section 5. In this 
case the related integration formula is [20] 


G(x,t) (l“t)'^(l+t)^dt r 2 Wj^G(x,t^) , 


-KRe(a,3)<l 


(6.44) 


where t, are the roots of 
k 


= 0 . k=l,..,n 


(6.45) 


and the weighting constants are given by 

( 2n+a+3+2 ) r (n+a+1 ) r (n+3+1) 2“'^^ 




(n+l):(n+a+3+I)r(n+a+3+l)P ^“’^^(t, ) 

n+1 k dt n k 


. (6.46) 


Analogous to the numerical integration methods developed in Section 5, 

the integr&l equation (6.41) and the condition (6.42) may now be expressed as 

n 
E 
k= 


^1 k j 


. .,n-l , 


? = A , 


(6.47a,b) 


where 


-KRe(a,3)<0 , P^^°^’^^(t^) = 0 , k= l,..,n , 


n 


(a+1,3+1), \ r\ *1 1 

P 1 (x.) =0 , 1 = l,..,n-l , 

n-1 3 ’ ’ ’ ’ 


(6.48) 


and are given by (6.46). Equations (6.47) provide n equations to deter- 
mi-ne g(t ^^) , . . ,g(t^) . 


If at one or both ends the solution is required to be bounded (hencej 
necessarily zero), (6.46) and (6.47a) will still be valid, the points tj^ 
describing the location and the number of the unknowns g(t^) will still be 
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obtained from (6.45) and the locations giving the number of equations 
will be obtained from 


(a) 0<Re(a)<l , 0<Re(3)<l : ° ^ j = l,..,n+l. 


(b) 

0<Re(ct)<l 

, -KRe(B)<0 


j = 1, . . ,n , 

(c) 

-KRe(a)<0 

, 0<Re(3)<l 

. p (a+l,B-l) = 0 , 

3 — l>**,n , 





(6.49) 

none of 

the cases 

given by (6.49) 

the condition (6.42) or 

(6.47b) is 


part of the formulation of the problem. In (6.49a) the additional equation 
(provided by , . . , ) is equivalent to the consistency condition of the 

integral equation and in (6.49b and c) the unique soluti<?n is obtained by 
simply solving (6.47a) v/ith j = l,..,n. As explained in detail in Section 5, 
the (physical) problem may, however, have additional conditions and unknown 
constants which may be handled in a straightforward manner. 

In the application of the numerical methods described in Section 5 and. 
in this section it is essential that special attention is paid to the con- 
vergence of the calculated results. At least one technique regarding the 
evaluation of the limit value (as n-»<») of certain calculated results is de- 
scribed in [25] and [27] . In the interest of space no numerical results 
will be presented in this article. However, there are certain techniques 
related to the methods described in this article V7hich are quite useful for 
extracting relevant physical information from the solution. Description of 
such techniques also falls outside the scope of this article. Among these 
one may mention the methods used to evaluate such physical quantities as 
the stress intensity factors, the strain energy release rate, and the crack 
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opening displacement directly in terms of the (calculated) density func- 
tions of the integral equations (see, for example, [28], [29], [31]-[34]). 
The detailed numerical results of the torsion problem discussed in this 
section may be found in [35] and [36]. In addition to [28, 29, 35, 

36], further applications of the singular integral equations with general- 
ized Cauchy kernels may be found in [37-39] . 

7. SINGULAR INTEGRAL EQUATIONS OF THE SECOND KIND 

Some relatively very simple mixed boundary value problems in mechanics 
give rise to singular integral equations which are of the second kind. For 
example consider the following basic formulas for the elastic half plane 

y<0 relating the surface tractions and the displacement derivatives: 


4y 3 


1 +K 3 x 


I 00 Jj. 

u(x,0) = yi (x,0) + — / a (t,0) — — , -oo<x<oo 

i yy\ J / 7T -“ Xy t-X 


4jj 3 


1 +K 3 x 


v(x,0) = ~ ya (x,0) + -^ /„ cr (t,0) -7^ , -oo<x<oo 

' ’ ' I jjy V 3 _00 yy ' J ^ t— X ' 


^ 0 (x,0) 

2y yy^ 


3 ... 1 “ 3 , dt 

- Y T" u(x,0) - — f v(t,0) — — , 

3x 7T 00 3t t-x 


-oo<x<oo 


^ a^y(x.O) = Y v(x,0) - ^ u(t,0) ^ 


-oo<x<oo (7.1a-d) 


where u, v are, respectively, the x, y- components of the displacement vector, 
(i,j = x,y) are the stress components and y = (k-1)/(k+ 1). From (7.1) 
it is clear that if the displacement vector is specified on part of the 
boundary y = 0 and the traction vector is specified on the remainder, (7.1a 
and b) or (7.1c and d) give the integral equations of the problem which are 
coupled and are of the second kind. It is also clear that instead of a half 
plane if one is dealing with an axisymmetric or plane problem for a domain 
bounded by a "smooth" surface the system of integral equations of the related 
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mixed boundary value problem will have a dominant part which is identical 
to that given by (7.1). 

As a very simple example consider the plane contact problem with 
friction shown in Figure 7 in which a rigid stamp of given profile is 
pressed against an elastic half plane. Assuming the coefficient of fric- 
tion p to be constant and defining 


ayy(t,0) = - p(t) , a^y(t,0) = - np(t) , a<t<b , 

v(x,0) = f(x) , (a<x<b) , 0) = PY (7.2) 


Equation (7.1b) gives the integral equation of the problem as follows: 
b 


03p(x) - “ / — dt = f(x) , a<x<b 

^ 7 T a t-x 


(7.3) 


Here the input function f (x) is, aside from a constant multiplier, the de- 
rivative of the function describing the stamp profile (i.e., f(x) = F'(x), 
where y = F(x) is the stamp profile). At the end points if the contact 
is smooth, the constants a and b defining the contact area are unknown. 

In this case the solution (or the input function) has to be such that the 
follot^ing consistency condition is satisfied: 


b 

/ [o)p(x) 
a 


1. f° P(t) 

7T a t-x 



dx 

w(x) 


- f (x) 
a w(x) 



(7.4) 


where w(x) is the fundamental function of the singular integral equation (7.3) 
Also, if the load P is specified the contact pressure must satisfy the follow- 
ing equilibrium condition 

'■ b . ... 

/ p(t)dt = P . (7.5) 

a 
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In solving the problem the conditions (7.4) and (7.5) account for the un- 
knowns a and b. 

For an elastic medium with a more complicated geometry the foregoing 
formulation basically remains the same with the only difference that the 
integral equation (7.3) now contains a Fredholm kernel which comes from 
the change in the geometry. For example, in the symmetric contact problem 
for an infinite elastic wedge of angle 26^ shown in Figure 0, using the 
standard Mellin transform technique the integral equation for the contact 
pressure may be obtained as 


1 ^ 

wp(r) - - / [ 
'' a 


1 

t-r 


+ k(r,t)]p(t)dt = f(r) 


a<r<b , 


(7.6) 


where 


p(r) = - ap,o(r,0^) , f(r) = Ug(r,0^) = F’(r) , 


k(r,t) 


'00^"»"o^ ’ 4u 
1 


2 

1 „ sin 0 -nsinO cos0 

1 TT o ' o o 


r log(t/r) t-r r 20^+sin20Q 

CO sinh20 y 


+ / ( 1 - 
o 


o' \ cospy j 

V )ti ^ dy 

D(y) ' ' r 


03 Tisin20 +cos20 -cosh20 y . ^ 

o ^ sinpy 


(1 + 


D(y) 


dy , 


(7.7) 


D(y) = sinh20^y + y siti20^ , p = log(t/r) 


(7.8) 


The problem is formulated in polar coordinates and F(r) describes the 
profile of rigid stamp. 


7 . 1 The Fundamental Function 

The dominant part of the singular integral equations of the second kind 
is of the following general form 


(7.9) 


A(J)(x) + I / dt = f (x) 

^ 7T a t-x o 


a<x<b 


where the hounded function may contain the part of the integral equation 
with the Fredholm kernel . Defining 


1 




F(z) = ^ dt , 

7T a t-z 


(7.10) 


and using the Plemelj formulas (3.61), (7.9) may be reduced to the following 
Riemann-Hilbert problem for the sectionally holomorphic functioii F(z) : 


f'^(x) - ~~ F~(x) = 2if^(x)/ (A+iB) 


(7.11) 


Considering the corresponding homogeneous equation 


x’Cx) = 0 , 


(7.12) 


the fundamental solution X(z) and the fundamental function w(x) of (7.9) 
may be obtained as [3] 


X(z) - (z-b)*^(z--a) ^ , w(x) " (b-x)^(x-a)^ , 


a = log (^^) + N , 3- 


— i— lop- + M 

27Ti ''A+iB^ 


2TTi ® 'A+iB 

where N and M are arbitrary (positive, zero, or negative) integers 
(7.11) and (7.12) we have 
+ 


(7.13) 

( 1 ) 


From 


/ F(x)v ' _ ^ F(x)\ _ 

^X(x)'^ 


the solution of which is 


2if (x) 

Q 

(A+iB)x’^(x) 


a<x<b 


(7.14) 


F(z) 

X(z) 


1 ^ 

TT a 


.f^(t)/(A+iB) 
(t-zlx"^ (t) 


dt + C 


(7.15) 


( 1 ) 


Note that if A and B are real then the exponents a and 3 are also real. 


For an example of complex exponents see Section 3,3. Also, see [4Q] for 
the equivalence of the two formulation. 
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where C is an arbitrary constant. From (7.10) and (7.15) the solution of 

(7.9) may be expressed as 

/ N b f (t)dt . 

Cj}(x) = - ^ + C„w(x) + ^5 f ^<x) , 


7T a (t-x)w(t) o 




a<x<b 


(7.16) 


where (7.13) is used to replace X (x) by w(x) and 


1-1 mi-i / / . . \ 

C = - BCe / (A-iB) 
o 


(7.17) 


The index of the integral equation is again defined by 


K = - (a+3) = - (N+M) 


(7.18) 


The general remarks made in Section 5 in connection with the contact problem 
shown in Figure 4 and the singular integral equation of the first kind re- 
garding the determination of the arbitrary integers N and M or the index of 
the problem and the constant which appears in the solution (7.16) are 
valid for the integral equations of the second kind also and will not be 
repeated here. 

It should be pointed out that the fundamental function of the problem 
can also be determined directly by applying the method described in Section 6.2 
For this we let the solution of (7.9) be 


<|)(x) = g(x)w(x) , w(x) = (b-x)°^(x-a)^ 


(7. 19) 


From (7.10), (7.19), (6.17) and (6.18) it is clear that (6.21) is still valid. 

Thus, substituting from (7.19) and (6.21) into (7.9) and multiplying both 

— B ~~0L 

sides first by (x-a) ^ and letting x-^-a and then by (b-x) and letting x^b 

we obtain 
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cotTTg = A/B 


cotTra = - A/B 


(7.20) 


From 


cotirG = cot7r(e+K) , K = 0,+l, 


- Arccot(A/B) , 


TT 


(7.21) 


it is seen that (7.13) and (7.20) are identical. 

To develop a mimerical method for the solution of the singular integral 
equations of the second kind, it will again be assumed that (f)(t) = g(t)w(t) . 
At a first glance the exact solution given by (7,16) appears to be not of 
this form. However, the integral in (7.16) can be evaluated and it can be 
shown that the expression for (j)(t) indeed has no regular terms. To do this 
first observe that 

b f (t)dt b f^(t)dt 

w(x) / -r- — ' = X (x) / 

a (t-x)w(t) a 


$(z) = / 


b f^(t)dt 


(t-z)X (t) i-x7x" 


(t-xyx"^ (t) 

b f (t)dt a f (t)dt 

If ^ j. _o j 

(t-z)X' (t) ^ (t-z)X (t) 


AJ.--C f (T)dT 

A+xB j- o ' 


2iB s (t-z)X(t) 


(7.22) 


where S is the contour shown in Figure 9 and (7.12) and (7.13) have been used. 
Now if fQ(T) is holomorphic outside S and continuous up to S, the contour 
integral may be evaluated as follows : 


- / 


f (T)dT 


fo(z) 


P(z) 


2iTi g (t-z)X(t) X(z) 
where P(z) is such that for large |zj 


(7.23) 


fo(z) 

X(2) 


= P(z) + 0(l/z) 


(7.24) 
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On the other hand using the Plemelj formula, from (7.22) it follows 


that 


b f (t)dt 
f 2 , 

^ (t-x)x'*'(t) 


j [$^(x) + $'(x)] 


(7.25) 


Thus, substituting from (7. 22)^(7 .25) into (7.16) and observing that 


(A+iB)x'^(x) = e^^“w(x)(A^+B^)^^^e“‘^^“ 


(7.26) 


we obtain 


(f)(x) = [(aW)“^^^P(x) + C^]w(x) . 


(7.27) 


7.2 Solution by Orthogonal Polynomials 

Consider the following integral equation in the normalized interval (-1,1) 


A(|)(x) + ^ dt + /J k(x,t)(})(t)dt = f (x) , -Kx<l (7. 


28) 


Assume that the fundamental function and the index of the problem have been 
determined and are given by (see Section 7.1) 


w(t) = (l-t)“(l+t)^ 


K = - (o+e) (Q,+l,-l) . 


(7.29) 


Observing that w(t) is the weight of Jacobi polynomials, we may express the 
solution of (7.28) as 


(j)(t) = g(t)w(t) , g(t) = S c P ^“’^^(t) 

o n n 


-l<t<l 


(7.30) 


where the coefficients c^ are unknown. Consider now the follcwing property 
of Jacobi polynomials [23, 41, 26]: 


AP ^“»®>(x)w(x) + - P ^®’^^(t) dt 

n 7T -1 n t~x 


2"*^ B ^ (-ot,-3) 

siiiTra n-K 


(x) , -1<X<1 , 





Re(a)>-1 , Re(g)>-1 , Re(a) (0, 1 , . . ) 


(7.31) 


From (7.30), (7.31), and (7.28) it follows that 




k^(x) = /j^ k(x, (t)w(t)dt . 


(7.32) 


An effective way of solving for c in (7.32) would be to expand both 


n 


f— fit — 

sides into scries of Jacobi polynomials (x) , and to compare the 

coefficients. Thus, using the orthogonality relations (5.72) we obtain 


p”'^T> 00 

— : — — 9.(-a,-3)c., + S d, c = F. , i = 0,1. 

sxn-fra j" ’ J+K inn j ’ 


(7.33) 


p/ ^\x)k^(x)dx/w(x) , 


n 


^^ (x)f(x)dx/w(x) 


(7.34) 


It is again v7orthwhlle to consider the following three cases separately. 


(a) K = 1 

In this case solving (7.33) by the method of reduction [24], if the 
series is truncated at n = N, the system of N+1 linear equation will con- 
tain the unknovTOs c^, ..., . The additional equation to solve the problem 
is provided by the foi.lowing condition: 


.1 


/j (|)(t)dt = 


(7.35) 


which, substituting from (7.30) and using the orthogonality conditions, may 

, , ( 1 ) 

be expressed as 


^^^Hote that for i< - 1 the arbitrary constant appearing in the exact solu- 
tion (7.16) or (7.27) is also determined from (7.35). 
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0^00 (o.g) = 


(7.36) 


(b) K = 0 

In this case, truncateci at n = N, (7.33) has N+1 equations and N+1 
unknowns and gives the unique solution. 


(c) K= - 1 

/ ct 6 ^ 

In this case observing that (t) = 1, the first equation (j = 0) 


in (7,33) may be expressed as 


( 1 ) 


°° 1 d-w 1 Hv °° 

F -Ed c = /, f (x) -%r - /, E c k (x) 

o Q on n -I ' w(x) -1 w(x) q P 


= / 


dx 


-1 w(x) “ -1 ^(x,t)^(t)dt] = 0 


(7.37) 


which is the consistency condition of the integral equation and is seen to 
be automatically satisfied if this technique is used. If the boundary value 
has a symmetry with respect to x in the sense that, for example, in Figure 7 
a = -b, then there is only one additional unknown, b and (7,33) and (7,35) 
provide N+2 equations to determine c^,,,,Cj^, and b. However, if the problem 
has no symmetry (which is always the case when the friction is involved), 
both a and b are unknown. Theoretically, the consistency condition with 
(7.35) provide the additional equations to determine a and b. In this case 
an extra equation may be gained by writing (7.33) as 


2B 


N-1 


-4^ e.(-a,-6)c; . + E d. c = F. , j - 0,..,N , (7.38) 

sinTra 3-I n=o 3^ ^ 3 

which, with (7.35) provide N+2 equations to determine c^, .. ,c^__j^, a, and b. 

It is seen that if k(x,t) is zero and f(x) is a polynomial of finite 
degree, by expressing it in terms of a series of Jacobi polynomials ^\x) 


( 1 ) 


Note that (7.32) starts with ^^(x) . Therefore in (7.33) c_^= 0. 

ioo 



and using (7.31) the solution of the integral equation (7.28) may be 
obtained in closed form by simple observation. Such solutions of some 
simple examples regarding the contact problem shown in Figure 7 may be 
found in [26] . The solution of the wedge problem shown in Figure 8 with 
and without a crack initiating at the apex will appear elsewhere. The 
application of the technique described in this section to singular integ- 
ral equations of the second kind with complex coefficients may be found 
in [31], [33], [42] and [43]. Extension of the technique to a system of 
equations in which A, B, and k(x,y) are square matrices is given in [44] 
and [25] . 

7.3 Solution by Gauss-Jacobi Integration Formulas 

In order to solve the singular integral equation of the second kind 
(7.28) in a direct way one needs an integration formula for the dominant 
part 

K((j)) = A4>(x) + I /J ^ dt . (7.39) 


Following closely the procedure outlined in Section 5 (see also [26]) such 
an integration formula can indeed be developed [45], Analogous to (5.53) 
first consider the expansion 


p( ^\x)pf“’^^(x) - P^“’^^(x)pf~'^’"^^(x) 
n-K 1 n q-K ^ ' 

p^a,3)(^) 


n 


= - S 


k=i 


(7.40) 


xjhere 


,(a,3) 


n 


(t^P =0 , k = 1, . . ,n , 


(7.41) 



rp( ff- ')p I / p^^J 3) 

'■ n-ic ^ V 3 dt n 



(7.42) 
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If we now select x as the roots of 
r 


= 0 » r = l,..,n-K , 


(7.43) 


from (7.40) it follows that 


p(-a,-3), V ^ 5 n-K 1 

^ k=l ^ .n, R^ 


Ap(a,3)(t ) 

dt n 


t, -X 

k r 


(7.44) 


We next consider (7.30) and assume that g(t) can be approximated to a 
sufficient degree of accuracy by 


g(t) 2 2 c P^°^’^^(t) 

j=0 J J 


(7.45) 


From (7.30), (7.31), (7.39), (7.44) and (7.45) for n>p it then follows that 




8in7ra Apfa.6)(t , 

dt n k 


(7.46) 


Again, the only approximation in (7.46) is due to the truncation in (7.45). 
Using the properties of Jacobi polynomials (7.46) can be put into the stand- 
ard form [45] 


n 


Bg(t^) 


KlKxpi . 2^^ — 




k \ 

2”'’'‘^(2n-^a+g)^(n+a)^(n-l-e) 

Trn;r(n+a+B-H)p'“i®ht^) £ 

2°^~^^'Hr(a+n+l)r(3-tn-ll) 
Fn:r(n-m+3+l)(l-tJ) p^'^v^^t^)]^ 


(7.47) 


• ( 2n+a4-R+2) V (g+n+l ) 2 ( B+n+1) 

F(n-I-L) ! r(n-l-cri-g+l)P^®’ *^^ (t^)-^ (t^^) 


^--(7.48) 
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Using (7.47) the integral equation (7.28) may now be expressed as 


n Bg(t^) 


^ — — + Trk(x ,t, )3 = f(x ) , r = 1,. .,n-K 

k=l r ^ 


(7.49) 


where t and x are given by (7,41) and (7.43). For K = 1, (7.49) gives n-1 

tC IT 

equations to determine g(t^^) , . . ,g(t^) . The nth equation is obtained from 
(7,35) which becomes 


n 


k=l 


(7.50) 


For K = 0 (7.49) contains n equations and n unknowns and gives a unique 
solution. In this case (7.35) may be used to determine the unknown end value 
(e.g., a or b in Figure 7). 

For K = -1 (7.49) and (7.35) give n+2 equations to determine 
g(t j^) , . . ,g(t^) , and the unknown end values a and b. It can again be shown 
that by using this method the consistency condition of the integral equation 


[f(x) - k(x,t)(|)(t)dt] 


dx 


dx 


/ \ = f-, K[(j>(x) / V 
w(x) -1 w(x) 


(7.51) 


is automatically satisfied. This can be seen by observing that l/w(x) is 
the weight of ^\x) , K= -I, P^ ^\x) = 1, 


n 


K[(|)(x)l = - 


2B 


P 

E 


(-a, 3) 

simra j j+1 


(7.52) 


and by using the orthogonality conditions for the Jacobi polynomials, 


As ‘ an example consider the wedge problem sho\m by Figure 8 and form- 
ulated by (7. 6-7.8). Assume that the rigid wedge has flat faces with sharp 
corners at the ends a and b. Thus, (7.6) must be solved for f(r) = 0 and 
under the condition that 


p(r)dr - P = P^/[2(sin9^ - qcos9^)] - A^ 


(7,53) 
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Let the coefficient of friction be n = 0.5. From (7.13) it then follows 
that 

a = - 0,545167 , 3 = - 0.454833 , K = 1 , 

p(r) = g(r)w(r) , w(r) = (b-r)“(r-a)^ . (7.54) 

The solution is obtained in a very straightforward manner by considering 
(7. 6-7. 8) with (7.49), (7.50), and (7.54). 

Of particular interest in problems of this type is the strength of 
stress singularity around the singular points a and b . This is known as 
the stress intensity factor and, in this case, may be defined by 

k(a) = lim /2 (r-a) ^p(r) , 

r->a 

k(b) = 1 iig (h-r) %(r) . (7.55) 

Some sample numberical results are shown in Table 7.1 where c = (b-a)/2 and 
P is defined by (7,53), P^ being the total wedging force. 
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